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Abstract 



Lorentz' objection to the Uhlenbeck-Goudsmit proposal of a spin magnetic moment 
of the electron, namely that the electron's equatorial rotation speed would exceed the 
speed of light by a factor ~ 10, mutated into an objection against Lorentz electrody- 

^ ■ namics (LED) itself when the spin magnetic moment became established. However, 

^ . Lorentz' renormalization calculation, based on the early 20**^ century's notion of a 

H ' purely electromagnetic electron, does not qualify as proper from a modern perspective. 

This paper shows that renormalization treated properly does lead to a mathematically 
consistent and physically viable LED. A new, relativistically covariant massive LED is 
presented in which the bare particle has a finite positive bare rest mass and moment 
of inertia. The particle's electromagnetic self-interaction renormalizes its mass and 
spin. Most crucially, the renormalized particle is a soliton: after any scattering pro- 
cess its rest mass and spin magnitude are dynamically restored to their pre-scattering 
values. This guarantees that "an electron remains an electron," poetically speaking. 
A renormalization flow study of the limit of vanishing bare rest mass is conducted 
for this model. This limit yields a purely electromagnetic classical field theory with 
ultra-violet cutoff at about the electron's Compton wavelength! The renormalized limit 
model matches the empirical electron data as orderly as one can hope for at the level 
of Lorentz theory. In particular, no superluminal equatorial speeds occur. 

(c) (2000) The authors. Reproduction for non-commercial purposes of this article, in its 
entirety and by any means, is permitted. 



1 Introduction 

The roots of the mass renormahzation program of quantum electrodynamics (QED) trace 
back to the late 19*'^ and early 20**^ centuries' notion that the electron's inertia has a purely 
electromagnetic origin. This notion was central to the pre-quantum attempts of the time to 
construct a consistent electromagnetic theory of matter and radiation that would unify elec- 
trodynamics and mechanics, now known as classical electron theory. When Lorentz proposed 
that Abraham's [Q, |^ purely electromagnetic, strictly rigid spherical electron (see Appendix 
A. 2) had to be replaced by an electron that undergoes the Lorentz-Fitzgerald contraction, the 
laws of special relativity began to take shape in the papers of Lorentz 0, and Poincare [Q; 
see [0 for an excellent account of the historical developments. Ironically, special relativity 
and renormahzation later contributed to the demise of relativistic Lorentz electrodynamics 
(LED) more than quantum mechanics. Of course, quantum mechanics removed LED from 
the list of contenders for a truly fundamental electromagnetic theory of matter and radia- 
tion, but it did not jeopardize the viability of LED as a physical theory in the classical limit. 
After all, Bohr's correspondence principle was undisputed in the quest for the fundamental 
quantum theory, and as such one would not have been, nor would one be unhappy with a 
consistent classical microscopic relativistic theory of electromagnetic matter and radiation 
that could serve as the empirically correct classical limit of a consistent quantum theory. 
Special relativity on the other hand had emancipated itself in the works of Einstein ||^, |^ 
and Minkowski P] as a theory of spacetime structure independent of LED or any other 
particular matter model [^. Thus, when Lorentz objected to the Uhlenbeck-Goudsmit [Q 
proposal of an electron spin magnetic moment by pointing out that LED would forbid the 
required magnitude of about one Bohr magneton because the purely electromagnetic electron 
would have to reach an impossible equatorial rotation speed of about 10 times the speed of 
light (see [jlO|, P- 35), he had furnished an argument that — eventually — weighed heavily 



against the physical sensibility of LED itself, and which has persisted ever since. 

Judged from our modern perspective on mass renormahzation ||TT1 , |12| , the purely electro- 
magnetic calculations by Lorentz ||^, |13|] and his contemporaries do, however, not qualify as 
proper. A proper renormahzation analysis would investigate the 'flow' toward vanishing bare 
rest mass for a consistent relativistic LED with positive bare rest mass ('massive' LED) in 
which the bare particle parameters of the model can be chosen such that the particle charge, 
the particle spin magnetic moment, and the renormalized particle mass match the empirical 
electron data without involving superluminal gyration speeds of the model electron. If the 
limit of vanishing bare rest mass exists, this limit would define a renormalized LED which in 
a sense completes Lorentz' program of a purely electromagnetic electron theory, and whose 
physical viability is no longer objectable to on basis of the superluminal arguments that 
grounded the original Lorentz theory. 

In this paper we will show that such a properly renormalized approach to the Lorentz 
program of electrodynamics does indeed lead to a satisfactory Lorentz electrodynamics. For 
simplicity we consider only a single Lorentz electron coupled to the Maxwell-Lorentz fields in 
the whole spacetime. The extension to a many particle theory does not pose any difficulties. 

We first present a novel, manifestly Lorentz-covariant massive LED which displays several 
crucial features expected of a realistic, consistent classical electrodynamics. In particular: 



• the dynamical equations constitute a Cauchy problem for the evolution of the physical 
state of massive LED; 

• scattering does not change the nature of the renormalized particle, which moves as a 
soliton; 

• the parameters of the bare particle can be chosen such that charge, magnetic moment, 
and renormalized mass of the 'free' particle can be matched to the physical electron 
data without involving superluminal gyration speeds. 

We then study the renormalization flow to vanishing bare rest mass of our empirically 
matched massive LED. The 'renormalized LED' which emerges in the limit is characterized 
by the following additional features: 

• the equatorial gyration speed reaches the speed of light c precisely in the limit of 
vanishing bare rest mass so that the bare gyrational mass converges to a nonvanishing 
contribution. It is fitting to call this a 'photonic' mass; 

• the renormalized spin magnitude is a derived property in our model. It turns out 
to be 3h/2 rather than h/2, with corrections of order a (Sommerfeld's fine structure 
constant), and we have g ~ 2/3 rather than g ^ gQ,^ = 2 or g ^ g^ = 1. 



• our renormalized LED is, strictly speaking, a classical field theory with ultraviolet 
cutoff corresponding to the physical electron's Compton length! 

It is not unexpected that in our classical theory the relation between spin and magnetic 
moment holds with a Lande factor g ^ 2. On the other hand, since we operate here with a 
minimalist model, it is not a priori clear whether a slightly more intelligently constructed 
classical model would not account also for the correct g ^ 2. In any event, the issue of the 
g factor should not be confused with the superluminal speed problem of Lorentz. 

Among the many physically correct features of the model, the truly rewarding one is its 
soliton dynamics. Without it the rest mass, magnetic moment and spin of the renormalized 
particle would typically be altered by each scattering process, and the theory would then 
simply fail to account for the classical motion of electrons and their electromagnetic fields. 

Another gratifying feature is the size of the renormalized, purely electromagnetic model 
electron, whose radius comes out to 1.5 physical electron Compton lengths, about 200 times 
larger than the somewhat misnamed "classical electron radius" e'^/mc^. Our renormalized 
LED is thus intrinsically equipped with the ultraviolet cutoff at which heuristic discus- 
sions [14 1 have traditionally set the spatial limit of applicability of classical electrodynamics 
as a physical theory. 

Remarkable is furthermore that the Cauchy problem for the state evolution bears a certain 
resemblance to the problem of motion of a black hole in general relativity. Unexpectedly, 
massive LED thus provides a playground for this much harder problem of general relativity. 



Our model owes much to the monumental work of Nodvik |jT5[ . However, although in Nod- 
vik's Lorentz-covariant massive LED one can match the particle parameters to the empirical 
electron data for charge, mass, and magnetic moment without encountering superluminal 



gyration speeds, one can do this only for sufficiently large bare rest mass. More specifically, 
the renormalization flow toward smaller bare rest mass terminates at about 1/2 of the em- 
pirical electron mass, i.e. before reaching the purely electromagnetic limit of vanishing bare 
rest mass, and its spin is then merely about 10~^h. In addition to this deficiency, Nodvik's 
massive LED suffers from a serious dynamical defect which it shares with the purely elec- 
tromagnetic models of Abraham and Lorentz; namely, the initial value problem is singular. 
Interestingly, both deficiencies of Nodvik's massive LED are caused by his choice of bare rest 
mass distribution, which concentrates all the bare inertia at the center of the particle so that 
its bare moment of inertia vanishes. By endowing the particle with a bare rest mass distri- 
bution that gives it a strictly positive bare moment of inertia our massive LED overcomes 
both deficiencies and, in this sense, improves decisively on Nodvik's massive LED. 

The presentation of our results is organized as follows. To begin with, in section 2 we 
stipulate the relativistic notation used in this paper. In section 3 we discuss the relativistic 
kinematical prerequisites of the model. Massive LED is defined by a Lorentz invariant action 
functional that 'minimally couples' the dynamics of the bare particle and the electromagnetic 
fields. We ffist formulate the relativistic mechanics of a gyrating bare particle for given 
Minkowski force and torque of unspecified origin (section 4) and then, in section 5, present 
Nodvik's []T3| covariant Maxwell-Lorentz equations which determine the electromagnetic 
fields also in our model, given the kinematical history of the charged particle. Our action 
functional for massive LED is introduced in section 6, where we also discuss its symmetries 
and associated conservation laws. The dynamical equations of our massive LED are presented 
in section 7. In section 8 we address the Cauchy problem for the state evolution. The 
scattering problem is discussed in section 9, where we will see that the renormalized particle 
behaves like a soliton. In section 10 the renormalization group flow toward vanishing bare 
rest mass is studied, in detail for the stationary particle, and in an outline for the dynamical 
particle. We conclude the main body of our paper in section 11 with a summary and an 
outlook on follow-up work and open problems. In particular, we also speculate about possible 
implications of our results in the quest for a consistent QED. 

In order not to interrupt the flow of the presentation, some technical material adapted 
from [1^, e.g. the derivations of the Euler-Lagrange equations, has been relegated to Ap- 



pendix A.l. To facilitate the comparison with Nodvik's model [|T5|, Appendix A. 2 contains 



a brief discussion of the point mass limit in which our massive LED formally reduces to 



the massive LED of Nodvik |15|. In particular, we explain that this limit is singular and 
severely overdetermines the initial value problem. Incidentally, a similar assessment holds 
for Abraham's |]l], ^ purely electromagnetic versions of semi-relativistic LED with and with- 
out spin. Since the ill-posedness of the initial value problem for the purely electromagnetic 
Abraham model is itself a new observation. Appendix A. 3 gives a concise summary of the 
semi-relativistic massive LED and its purely electromagnetic, singular limit. There we start 
with the massive model with spin and show next that the more familiar semi-relativistic 
massive 'LED without spin' actually is a theory with spin in 'hidden' form, obtained from 
the massive model with spin in a non-singular limit. We then address the singular purely 
electromagnetic Abraham limit. A fair list of references has been selected from the vast 
pertinent literature. 



2 Minkowski space: vectors, tensors, and all that 

Various equivalent conventions for the Minkowski space formalism exist in the literature. By 
and large, we here follow the conventions of the book by Misner, Thorne and Wheeler [|T6|. 
However, since the model is essentially a flat space model, we will not introduce forms. 
Instead, throughout this work we use four-vectors and tensors constructed from them. 

Thus, abstract Minkowski space M^ is a fiat four- dimensional pseudo-Riemannian man- 
ifold with Lorentzian metric of signature +2. The invariance group acting on Minkowski 
space is the inhomogeneous Poincare group; its homogeneous subgroup is called the Lorentz 



group; for good discussions, see [|T7[ |l^. We can choose any point P G M^ and identify M^ 
with the tangent space Tp(M^) at P. A tetrad JFl = {eg, ei, 62, 63} of fixed unit four-vectors 
in Tp(M^) that satisfy the elementary inner product rules 

— 1 for yU = z/ = 
e^ ■ e,, = ■^ 1 for /i = z/ > (2.1) 

for /i 7^ z/ 

is a basis for Tp(M^), called a Lorentz frame. With respect to a Lorentz frame, any real- 
valued vector V G Tp(M^) has the representation 

3 
v = J2^^e, (2.2) 

At=0 

with f '^ G M given by f ° = — v ■ Bq and by f '^ = v • e^ for fi = 1, 2, 3. Therefore, Tp(M^) can 
be identified with M}'^, the set of ordered real 4-tuples v = {v^, f ^, v"^, v^), which is a global 
coordinate chart for M^. 

To emphasize spacetime among the various physical realizations of abstract Minkowski 
space, spacetime points will be denoted by E and called events, and their four-vector repre- 
sentation w.r.t. JFl will be denoted as x or y etc. The coordinates of x in a Lorentz frame 
JFl are sometimes grouped as 

X = (ct, a:;) , (2.3) 

where x = {x^,x'^,x^) is called a 'point in space,' and t = x^/c an 'instant of time,' with 
respect to the given Lorentz frame. Here, c is the speed of light in vacuo. Event four- vectors 
X are classified into spacelike, lightlike, and timelike according as x • x > 0, x ■ x = 0, 
or X ■ X < 0, respectively, where the inner product of any two event four- vectors x and 
y is evaluated using ( p.2|) , distributivity and the elementary inner product rules ( p.l|) . A 
spacelike four- vector is connected through an orbit of the Lorentz group with a four- vector 
of the form (0, x), a timelike one with a four- vector of the form (ct, 0). The lightlike vectors 
form the 'light cone,' a double cone in M^'^ which separates the spacelike from the timelike 
four-vectors. This classification is then carried over to four-vectors v in abstract Minkowski 
space. 

For convenience, v ■ v will sometimes be abbreviated thus, 

v-v=\\vf (2.4) 



where || v|| is the Minkowski 'norm' of v, defined as the principal value of (v ■ v)^/^. Notice 
that ||v|| is negative for timelike vectors. 

The tensor product between any two four- vectors a and b is a tensor of rank two, denoted 
by a (g) b, and defined by its inner-product action on four-vectors, thus 

{a^b)-c = a{b-c) (2.5) 

c-{a(g)b) = {a-c)b (2.6) 



(Notice that our convention here differs from the one used, for instance, in [IS].) Given a 
frame {e^}, any tensor of rank two, T, can be uniquely written as 

T= ^ r-e^®e,. (2.7) 

A rank-two tensor T for which T^'^ = ±T^^ in ( |2.7| ) is called symmetric {+ sign), respectively 
anti- symmetric (— sign). For symmetric (+) and anti-symmetric (— ) tensors we have 

T±-v = ±v-T±, (2.8) 

but for a general tensor 

' general ' *^ 7^ It V • I general • I'^'^J 

A special example for an anti-symmetric tensor of rank two is the exterior product be- 
tween two four-vectors, denoted by a wedge-up product a A b, and defined by 

aAb = a®b-b®a. (2.10) 

Examples for symmetric tensors of rank two are, (i) the symmetrized tensor product between 
two four- vectors, denoted by a wedge-down product aV b, and defined by 

ay b = a(^b+b®a, (2.11) 

and (ii) the metric tensor 

S= Yl a^'e^^e, (2.12) 

with g'^'^ = e^- Bij given in ( |2.1| ). Notice that g has the same components g^'^ in all Lorentz 
frames because the e^ ■ e^ are Lorentz scalars. Notice in particular also that g acts as identity 
on four-vectors, i.e. g ■ v = v. 

The (anti-) commutator of any two tensors of rank two, A and B is defined as usual, 

[A,B]^ = A-B±B-A. (2.13) 

The action of tensors on a vector is associative, i.e. 

(A-B)- v = A-(B- v) (2.14) 
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so that 

[A, B]^- v = A-(B- v)±B-(A- v). (2.15) 

Warning: In general, 

A-(vB) ^ (A- v)-B. (2.16) 

The four-trace, or contraction, of a tensor of rank two is given by 

TrT= ^ g^'^T^''. (2.17) 

In particular, Trg = 4, while Tr A = for any anti-symmetric tensor A. 

For a differentiable function /(x) we denote its four-gradient with respect to the metric 
g by Vg/. In local Lorentz frame coordinates it is given by 

VJ(x) = (-9J,V/). (2.18) 

where V is the usual three-gradient. In analogy with the conventional curl, we also define the 
four-curl of a differentiable four- vector function as the anti-symmetric four tensor function 

VgAA(x) = Yl ^'''^''e^®e,(eA-Vg)(e^-A) (2.19) 

where the e'^'^^^ are the entries of the conventional rank-four Levi-Civita tensor. The four- 
Laplacian with respect to g is just the (negative of the) d'Alembertian, or wave operator, 
given by 

Ag = Vg ■ Vg = -D (2.20) 

3 Particle kinematics 

A spherically extended, 'rigidly' gyrating particle in spacetime is kinematically described by 
a pair of curves in Minkowski space. We shall assume that both curves are of class C^. 

One curve is the so-called world-line of the particle, a timelike future-oriented curve in 
spacetime itself. The concept of the world-line is familiar from the relativistic mechanics of 
point particles. For an extended 'rigid' particle an event on its world-line marks the location 
of the particle's center in spacetime. The world-line has associated with it an invariant 
'arc-length' element dr. In units such that the speed of light c = 1, which henceforth we 
shall use, this invariant arc-length element becomes the proper-time element. Choosing any 
particular event Eq on the world-line as reference world point and assigning to it the proper- 
time t{Eq) = 0, and stipulating the 'future' direction as the one along which dr increases, 
then by integrating dr along the world line automatically determines a unique invariant 
proper-time t{E) G M for any other event E on the world-line. Since the world-line is 
timelike future-oriented, its proper-time assignment is one-to-one onto and order-preserving. 



whence serves as a natural parameterization for the world-hne. In any particular Lorentz 
frame JFl, the particle world-line is then given by a mapping r i— > x = q(T) G M^'^. The 
proper-time element accordingly is given by dr = -y/— dx ■ dx, where dx is taken along the 
world-line. 

With this parameterization of the world-line, the four- velocity of the particle at the event 
q(r) is obtained as dx/dr = q(r) = u{t). The map r i-^ u(r), called 'four-hodograph,' 
traces out a complicated curve in four-velocity space, located on the future unit hyperboloid 
which is the genus- zero three-dimensional sub-manifold defined by u ■ u = — 1 and u^ > 0. 
When r I— *■ u is constant, the four-hodograph degenerates into a Minkowski point. We will 
sometimes group the four-velocity into its time and space components w.r.t. the chosen 
Lorentz frame, as 

u = (7,7i;), (3.1) 

where v is the usual three-velocity in the space part of the chosen Lorentz frame, and where 

7= r^ (3.2) 

Finally, by taking the second derivative w.r.t. r we obtain the four-acceleration of the 
particle at the event represented by ^(t), thus d^x/dr^ = u{t) = b{t). The map r i-^ a{T) 
is a section in the tangent bundle to the future unit hyperboloid of four- velocity space. 

The world-line does not contain any kinematical information about the gyration of the 
extended rigid particle. That information is contained in what may be called the gyrograph 
of the particle, a spacelike-oriented curve in four angular velocity space, or gyrospace for 
short. The introduction of gyrospace and the gyrograph map requires the consideration of 
two non-inertial reference frames: 

• the body frame J?-body; 

• the Fermi- Walker frame J^yw- 

The body frame J^body is a co-moving and co-rotating frame, fixed in the particle. It is 
given by a tetrad of orthonormal basis vectors {e^}/i=o,...,3, with Gq = u(t), that satisfy the 
equation ||16[ 



— e^ = -ri-e^, (3.3) 

where Vt is the anti-symmetric tensor of the instantaneous rate of four-gyration of the par- 
ticle's body frame J^hoAj as seen in the Lorentz frame JFl. The components of four-vectors 
with respect to J-body will be marked by primes: x = X]o<u<3^^^m- 

The Fermi- Walker frame jFp„ is the natural co- moving tetrad of orthonormal basis vectors 
{e^}^=o,...,3; with eo = ij(t), satisfying the law of Fermi-Walker transport |T^, 



d 
dr 



-j-e^ = -fipw-e^, (3.4) 



where 

fipw = u A u (3.5) 

is the anti-symmetric Fermi- Walker tensor. The components of vectors with respect to jFp^ 
will be marked by overbars, i.e. x = X]o<u<3-^^^m- 

Fermi- Walker transport is a kinematical effect which does not discriminate between ob- 
jects of different physical characteristics like inertia or charge. The most prominent effect 



associated with Fermi- Walker transport is the celebrated Thomas precession |20| of the 



space axes of JFp^ relative to the Lorentz frame JFl, described by the angular velocity three- 



vector 10, EO, 21 



<^T = (7 — 1) 2^ ^ '^ • (^■6) 

\v\ 

Here a is the usual three-acceleration of the origin of jFp^ as seen from JFl, and v its three- 
velocity. 

If the four-gyration of the body frame w.r.t. the Lorentz frame differs from the Thomas 
precession, one ascribes this difference to the particle's intrinsic four- gyration, which is 
captured by the anti-symmetric tensor 

o - ripw = n^ . (3.7) 

The subscript '5;' here is meant as reference to Euler's monumental analysis of Newtonian 
gyration of an inert rigid body. Beside being anti-symmetric, the tensor fig satisfies 

n^u=0, (3.8) 

so that fig has only three independent components. It follows that Op is dual to a spacelike 
four-vector iVp which satisfies 

Op ■ wp = (3.9) 

and Wp ■ u = 0. In J^yw, Wp is given in the form Wp = (0, ojp), where ojp is the usual angular 
velocity three- vector which is directed along the instantaneous axis of body gyration in the 
spacehke three-slice of the frame ^pw. 

At each instant of proper-time r, the instantaneous state of intrinsic four-gyration of the 
particle is captured by the four- vector Wp(r). As r G M varies along the particle's world-line, 
the map r i— > We{t) ^ I^^'^ traces out the gyrograph of the particle in gyrospace, which is 
the Minkowski space consisting of all possible four-vectors with the meaning of an Wp. The 
gyrograph is spacelike, since for ojp 7^ 0, we have 

-iTr(rjp-Op) = llwpf = |a;p|2>0. (3.10) 

4 Bare particle mechanics 

Formally stripping the particle from its electromagnetic features leaves one with the bare 
particle. Our bare particle is characterized by a strictly positive inertia against changes of 
its translational and gyrational momenta, respectively. 
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4.1 Bare inertia 

In our model, all forms of bare inertia are computed with a compactly supported, spherical 
'density' /m(| . |) : M.^ -^ M+ of inert mass of unspecified non-electromagnetic origin. In 
general, /m : IR^ -^ 1R+ is a measure on [0,R]. We demand that /m generates a strictly 
positive moment of inertia (see below). For instance, two admissible mass densities are: 

3 

• (volume inertia) /m(||x||) = mt— — Xfl(||x||), 

• (surface ineHia) /m(||x||) = m^——^5{\\x\\ - R) . 

Here, (5( . ) is the Dirac distribution, and Xr{ ■ ) the characteristic function of the inter- 
val [0, B\ with i? > 0. A fraction of the mass in /m is allowed to be a point mass, but not 
all of it. Notice that the complete point mass density /m(|a;|) = mb(l/47r|a;p)(5(|a;|) violates 
our condition that /m must have a strictly positive moment of inertia. The singular limiting 
case of a point mass density is addressed in Appendix A. 2. 
The most basic inertia is the hare rest mass, given by 

m^= [ U\x\)d'x, (4.1) 

where r = |a3|. Since in general the bare particle is in a state of gyration, its instantaneous 
effective inert mass in the Fermi-Walker frame is the sum of its bare rest mass and the 
relativistic kinetic energy associated with its Euler gyration, which we call the bare gyrational 
mass, a Lorentz scalar given by 

^b(l|wK||)= /" , ^ /^(||x||)^(u-x)d^x (4.2) 

l-lia-xf 



E 

R 



artanh(|a;E|r)/m(r)47rr"dr . (4.3) 



Kfp r 



The 'integration over a simultaneity space slice' in ( |4.2| ) follows the insight of Fermi [^, who 



used this setup in his solution to the infamous 4/3 problem (see ||23|, |2^, |2^, |^, ^, ^ ^ ^| 



for discussions of the 4/3 problem). In the limit of small inertial angular velocity, \lJe\ ^ 1, 
a MacLaurin expansion of ( [4.3|) in powers of \ujj^\ gives the expected result 



A^b(I^El) =m^ + ^/blc^Ep + 0(|a;,|^) , (4.4) 

where 

h = l [ \x\'U\x\)d'x (4.5) 

•J Jm.3 

is the bare particle's non-relativistic principal moment of inertia. By the spherical symmetry 
of the particle in any of its rest frames the principal inertias are of course degenerate. The 
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corresponding non-relativistic moment-of-inertia three-tensor of the bare particle generahzes 
to the tensor of the Minkowski inertia (about q), given by 

lb(||w,||)=/ ll^ll'g-^^^ ./^(||x||)^(u(r)-x)d^x. (4.6) 

'''' ■'i-mjT)-xf 



which captures the inertia associated with changes of the state of intrinsic four-gyration of 
the bare particle. 

4.2 Bare dynamics 

The various inertias acquire their meaning from their role as 'coefficients of resistance' against 
changes of the bare particle's various Minkowski momenta. The bare particle's four-vector 
of Minkowski momentum is given by 

Pb = 7Wb", (4.7) 

and the four- vector of its Minkowski spin by 

Sb = lb ■ We . (4.8) 

For a given Minkowski force f{T) and Minkowski torque t(r), here of unspecified origin, but 
satisfying t- u = and f ■ u = 0, the evolution of these momenta is governed by the coupled 
dynamical equations 

d „ , , 

^Pb = f (4.9) 

and 

d 

— Sb + f2Fw-Sb = t. (4.10) 

dr 

Due to the occurrence of q(r) in the gyrational mass, these dynamical equations are of 
second order for q and of first order for Wj^, so they need to be supplemented with initial 
data for q, u, and w,^. A discussion of these equations is interestingly in itself, but will not 
be pursued here. 

5 Nodvik's Maxwell— Lorentz equations 

If the charged particle's world-line ^(t) and gyrograph fl^iT) are given, the anti-symmetric 
electromagnetic tensor field of rank two F can be computed by solving the linear, relativistic 
Maxwell-Lorentz equations 

Vg-*F = 0, (5.1) 

Vg-F = 47rJ, (5.2) 
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where *F is the (left) Hodge dual of F, and where the covariant charge- current density four- 
vector field J is of the general form 'charge density times four- velocity.' The specific form 
of the charge- current density of course depends on the particle model under consideration. 



Our massive LED inherits the charge-current density from Nodvik's [0] model of a gyrating 
Lorentz electron, which generalizes the non-relativistic expression for a gyrating, strictly 
rigid Abraham electron (Appendix A. 3) to the manifestly covariant expression 

J(x) = T" I7(x; r) /e(||x - q(r)|| ) 6(u{r) ■ (x - q(r))) dr . (5.3) 

Here, 

l/(x; r) = u{t) - n{T) ■ (x - q(r)) , (5.4) 

is the four- velocity of a particle element located at x in JFl, with ft = Opw + ^^e, and 
/e(| . I) : M^ — > M~ is the spherically symmetric electrical charge 'density' (in general a 
finite measure) of the particle at rest. To avoid infinite electromagnetic energy densities 
and forces, we request that fe does not feature any singular Dirac point charge distribution 
anywhere on its support; more general Dirac distributions are however allowed. For instance, 
the following two possibilities are admissible: 

3 

• (volume charge) /c(||x||) = -e^^^X«(||x||), 

• (surface charge) /e(||x||) = -e— — 5(||x|| - R). 

Here, e > is the elementary charge as appropriate for a model of an electron. 

Problems with infinite self energies and self forces aside, however, it should be noticed 
that in the formal limit /e(|cc|) — *■ — e(47r|a;p)~^(5(|a:;|), together with the assumption that 
II We II < C, the Nodvik charge-current density (|0|) reduces to the familiar expression for 



the spinless point charge |]26|, ^ 



/ + 00 
u(r)5(x- q(r))dr. (5.5) 

-oo 

6 Minimal coupling of particle and fields 

In order to obtain the dynamical equations of massive LED, we have to couple the relativistic 
Maxwell-Lorentz equations for the electromagnetic fields and the relativistic equations for 
the particle's world-line and gyrograph in a consistent manner. Since both the dynamics 
of the bare particle and the dynamics of the Maxwell-Lorentz field derive separately from 
some principle of 'least' action, we can follow Nodvik ||T5[ and apply the standard proce- 
dure of minimal coupling to the action functionals and then postulate the principle of least 
action for this new, joint action. As a general reference for relativistic action principles we 



recommend |32 . 
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Figure 1: An example of an admissible four-domain S, which contains an ordered history of 
the particle. The extra-bold line is the world-line of the particle. The solid lines indicate the 
support of /e(||x — q(r)ll) restricted to the simultaneity space shces u{t) ■ (x — q(r)) = at 
Ti, T2 and an intermediate value of r. The dashed lines are the space slices that bound S. 




Figure 2: An example of a non-admissible four-domain S. No ordered particle history is 
possible because the supports of the functions x i-^ /c(||x — q(T)||) restricted to different 
space slices of simultaneity u{t) ■ (x — q(r)) =0 must intersect. 



6.1 The action of massive LED 

The action is defined through an integral over any four-domain S which is sandwiched 
between two disjoint spacelike three-subspaces that contain the intersections of the supports 
of /e(||x — q(T) II) (and of /m(||x — q(T)||)) with the spacelike slices u{t) ■ (x — q(T)) = of 
events which are simultaneous to the events at q(r) for r = ri and t = T2, respectively. Such 
a four-domain is called 'admissible.' It contains ordered particle histories, defined by particle 
world-lines which connect the supports of /e(||x — q(T)||) in the slices u{t) ■ (x — q(r)) = 
at r = Ti and r = r2 in such a way that successive space slices of simultaneity do not 
intersect in the support of /e(||x — q(r)||); see Figure |l|. Not every world-line corresponds to 
an ordered history. The space slices of simultaneity for two different values of r may intersect 
in the support /e(||x — q(r) ||) if the particle acceleration is too strong; see Figure |[ 

To formulate our minimally coupled action functional, we also need to introduce the four- 
vector of the electromagnetic potential. We recall that (p.lD implies that F is the exterior 
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derivative of an electromagnetic potential four- vector A(x), i.e. F(x) = Vg A A(x). We 
will work within the class of A that satisfy the so-called Lorentz gauge condition Vg ■ A = 0, 
which still allows the gauge transformation A ^ A + VgZ/' for any ip satisfying the wave 
equation Dip = 0. 

Given an admissible S, the action functional A is given as integral 



A= /:(x)d*x, (6.1) 

where C{x) is the Lagrangian four-density for the minimally coupled field-particle system, 

''-' •'-field ~r ^-coupl ~r •'-bare • \^"^) 

The first term is Schwarzschild's Lagrangian density for the free electromagnetic field of 
Maxwell, i.e. 

Aeid(x) = --^ Tr (F(x) . F(x)) , (6.3) 

and the second term is the standard Lagrangian density for the minimal electromagnetic 
coupling of particle and field degrees of freedom, i.e. 

/:coupi(x) = J(x)-A(x), (6.4) 



where J is the Nodvik charge-current density given in ( ^.3|) . The last term is the Lagrangian 
density for our bare particle, 

Cb.re{x) '=' - r ^l-||f2,-(x-q)||' U{ ||x - q|| ) 5(u ■ (x - q)) dr . (6.5) 

6.2 Invariants for the action functional 

Our action functional is invariant under several symmetry transformations of base and tar- 
get spaces, associated with which are invariants that give the physical conservation laws. 
These invariants are listed below. We consider only ordered histories, so that the domain of 
integration in the conservation laws can be written M^'^ instead of S. 

6.2.1 Invariance under base space transformation 

Invariance of the action functional under translations of spacetime implies conservation of 
the total Minkowski momentum four- vector P of the field-particle system, the constant time 
component of which expresses energy conservation, the space components the conservation 
of linear three-momentum. The Minkowski momentum of the field-particle system can be 
written in the form 

P = Mfp(r) • u{r) , (6.6) 



14 



where 



Mfp(r)=A<b(||w,(r)||)g+M,(r)+ / T(x) 4u(r) ■ (x - q(r))) d^x , (6.7) 

is the symmetric field-particle tensor mass, itself composed of the following tensor masses: 
the first term on its r.h.s. is just the bare gyrational mass of the particle written as diagonal 
tensor with the help of the metric tensor g; the second term on its r.h.s is the symmetric 

Nodvik tensor mass, 



M.(r) = - / [x X, [F(q(r) + x), f7,(r)] J f,{ ||x|| ) 6{u{t) ■ x) d^x ^ 

JlRi.3 



(6.8) 



extracted (see Appendix A.l) from the Minkowski momentum four-vector of electromagnetic 
spin-orbit coupling given in |T^; the last term on its r.h.s. finally is the electromagnetic field 
mass tensor, with 



T(x) = 1- f F(x) . F(x) - 1 Tr(F(x) ■ F(x))g) (6.9) 



An 



the conventional traceless tensor of the electromagnetic energy-momentum-stress density, a 
symmetric tensor field of rank two. It's four-divergence, given by 

Vg ■ T(x) = F(x) • J(x) , (6.10) 

does in general not vanish on the support of the particle densities /m and fe] in a nutshell, 
this is the origin for the need of the Fermi prescription of integration over space slices of 
simultaneity (cf. |2^, |27[|). 

Invariance of the action functional under spacetime rotations implies conservation of the 
anti-symmetric tensor of total Minkowski angular momentum L of the field-particle system, 

L = Sb(r) + q(r) AP+ [ xAfT(x + q(r)) ■ u{t)) 6{u{r) ■ x) d^x , (6.11) 

where P and T are given above, and where 

Sb(r) = - / X A —^M=^^U{\\x\\ ) 6{u{t) ■ x) d^x (6.12) 



is the anti-symmetric tensor of bare Minkowski spin (about q) of the particle. The tensor Sb 
is dual to the bare spin four- vector (^.8|), hence it is of space-space type w.r.t. u, viz. 



Sb-u = 0. (6.13) 

Since L is anti-symmetric, its constancy expresses six individual physical conservation laws. 
The constancy of the three space-space components expresses the conservation of the total 
angular momentum, and the constancy of the three time-space components correspond to 
the statement that the center of mass of the field-particle system moves with constant three- 
velocity along a straight line in the three-space slice of the Lorentz frame. 
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6.2.2 Invariance under target space transformation 

Invariance of the action functional under gauge transformations A ^ A + Vg^ implies 
charge conservation locally, i.e. (|5.iJ|) satisfies the continuity equation 



Vg-J(x) = 0, (6.14) 

which is demanded by the Maxwell-Lorentz equations and readily verified to hold also by 
direct computation with ( |5.3|) , cf. |15]. Of course, ( |6.14| ) implies the conservation of the total 



charge 

Q = - [ eo • J(x) (5(eo ■ X - to) d^x , (6.15) 

i.e. Q = —e in our model of the electron, for all to- 

In addition, there is a symmetry associated with the spherical rest symmetry of the 
particle, giving invariance of the action functional under shift of the Euler angle around the 
instantaneous axis of gyration. (For the definition of the Euler angles, see Appendix A.l.) 
This corresponds to the conservation of 

s' = -^Tt {S\t) . S\t)) , (6.16) 

where 

S^ = S+ [u® u,S]^ (6.17) 

is the space-space part w.r.t. u of the anti-symmetric tensor of total Minkowski spin, 

S(r)^='SKr) + Sb(r), (6.18) 

with 

Sf(r) '^ / (x - q(r)) A A(x)/e( ||x - q(r)|| ) dUr) • (x - q(r))) d^x (6.19) 

JRI.3 ^ ^ 

the tensor of the particle's electromagnetic Minkowski spin about q(T). 

We remark that despite its appearance, Sf is invariant under a gauge transformation 
A^ A + Vgip. 

The law for s"^ means that the dual four-vector to S evolves on the spacelike genus-one 
hyperboloid. 

6.3 The principle of 'least' action 

Our action principle is formulated for a fixed four-domain S, with the constraints 6u{t) = 
for r = Ti and r = T2. Likewise, the variations of the world-line Sq{T) and of the Euler angles 
(see Appendix A.l) vanish at r = ri and r = T2. Following standard physics terminology we 
refer to this Hamiltonian principle of stationary action as 'principle of least action,' although 
we do not investigate the second variation here. 
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The particle world-line satisfies the variational principle 

5(Aarc + Aoupl) = (6.20) 

under an independent deformation of the world-line r i— * q(T) (with r G [ti, T2]) into a world- 
line T \-^ q(T) + Sq{T) (with r in the same range), with Euler variables 6 and four-potential 
A fixed. 

The gyrograph satisfies the variational principle 

-^(Aare + Aoupl) = (6.21) 

under an independent variation of the gyrograph of the particle, with fixed world-line t ^—>- q 
and four-potential A. 

The fields satisfy the standard Schwarzschild principle |^ 0, 



5(Aeld + Aoupl) = (6.22) 

for prescribed world-line r i-* q and gyrograph r i-^ fig. 

7 The dynamical equations of massive LED 

The dynamical equations are the Euler-Lagrange equations for the variational principle 

6A = 0, (7.1) 

with q, A and the Euler angles treated as independent variables. We will state the rela- 
tivistic dynamical equations in a format closely reminiscent of the semi-relativistic equations 
of the massive Abraham theory described in Appendix A. 3. 

7.1 Field equations 

Recalling that F = Vg A A by definition, which already implies the homogeneous Maxwell- 
Lorentz equations ( p.l| ), the principle of least action together with Vg ■ A = now yields 
the equivalent potential form of the inhomogeneous Maxwell-Lorentz equations, i.e. the 
inhomogeneous wave equation 

nA(x) = 47rJ(x) , (7.2) 

with J Nodvik's four- vector of charge-current density given in (|5.3|) . 

7.2 World-line equations 

The world-line equations read 

^P = f, (7.3) 

17 



where 

p{r) = M(r) ■ u{t) (7.4) 

is the Minkowski momentum of the particle, with M = Ai^S + M^ (see ( |6.7| )), and where 

f{r)= f F(x)-U(x;r)/e(||x-q(r)||)<5(u(r).(x-q(r)))d^x, (7.5) 

is Nodvik's [|T5| Abraham-Lorentz type Minkowski force. 



7.3 Gyrograph equations 

The equations for the gyrograph read 

— Sb+ [ripw,Sb]_=t, (7.6) 

where 

t(r) = / (x - q{r)) A (F(x) ■ [/(x, r)) Ve( ||x - q(r)|| ) sUr) ■ (x - q(r))) d^x 

(7.7) 

is the anti-symmetric tensor of the Abraham-Lorentz type Minkowski torque, with a^ = 
(g + u ® u) ■ a the projection of the four-vector a onto the space shce w.r.t. u. The 
space projector g + u ® u under the integral guarantees that (|7.6|) preserves the intrinsic 
space-space character of Sb. The commutator with the Fermi- Walker tensor accounts for 
the four-precession of the Fermi- Walker-transported, natural co-moving frame with respect 
to which the inertial gyrations are defined. 

8 Evolution of the state in massive LED 

In the concise component-free Minkowski space notation the similarity of the world-line 
equation (|773|), supplemented with (|7.4| ) and (|775|), to the Newtonian type evolution equa- 
tion ( |A.84D of the semi-relativistic massive LED with spin is clearly recognizable, and so is 
the similarity of the gyrograph equation ( [7.6|) , supplemented with (fTT]), to the Eulerian evo- 
lution equation (|A.86|) . The almost familiar appearance of the dynamical equations of our 



relativistic massive LED is, however, deceptive. As we will now see, the physically natural 
interpretation of these dynamical equations as a set of evolution equations for the dynamical 
state of the theory unveils some radically new elements of physical reality (as defined by the 
theory) . 

The characteristic of a physical state, loosely speaking, is that its specification at any 
particular instant completely determines, via first-order evolution equations, the time deriva- 
tive of the state at that instant and, hence, allows to continue uniquely to the state at the 
next instant. To find out what precisely constitutes a 'state at any particular instant' in 



our massive LED, we have to analyze for which dynamical variables the equations pose a 
first-order initial value problem (Cauchy problem) in an orderly manner. This is done in 
the following two subsections. We then address the problem of nearly free global evolutions. 
Since the evolution equations of our massive LED are implicit in a rather unfamiliar manner, 
we first rewrite them into a quasi-explicit format. We then propose a direct generalization 
of the familiar Picard iteration to solve the equations — in principle. We finally remark on 
the question of global well-posedness of the Cauchy problem for arbitrary admissible states. 
This issue is only partially settled. Various possible 'mechanisms' that may cause some sort 
of singularity to occur at a finite time have to be controlled. However, for the special case of 
purely gyrational particle dynamics coupled to the Maxwell-Lorentz fields we succeeded in 
proving, by using our generalized Picard iteration, a global existence and uniqueness result. 



This proof is going to appear in 33 



8.1 The states in massive LED 

We first need to clarify the phrase at an instant of time. In special relativity, this is a quite 
flexible notion. Generally it refers to the collection of all the events on a three-dimensional, 
simply connected, spacelike hypersurface. We assume without loss of generality that such a 
hypersurface is defined as the level set S^ = {x : T(x) = t} of a differentiable function T on 
spacetime. Here, a three-dimensional hypersurface S^ is called spacelike if VgT(x) is timelike 
at each x in Sj. If the intersection of any two such spacelike hypersurfaces S^^ and S^j, with 
ti 7^ t2, is empty, and if furthermore the four- dimensional set UtgRSt covers spacetime, UtguSt 
is called a foliation of spacetime (generated by T), and each three-dimensional hypersurface 
Ej is called a leaf of that foliation. For example, the standard foliation of a Lorentz frame 
•^L = {^0, ei, 62, 63} into affine space shces of the form {x : a;° = t} is generated by the 
function Tjr^{x) = —Bq ■ x, which has a constant timelike four-gradient, ^gTyr^{x) = —eg. 
There are uncountably many foliations of spacetime among which one has to find any one 
that is compatible with the structure of the dynamical equations, preferably a convenient 
one. Unfortunately, as we shall see in a moment, the standard foliation is not compatible 
with the dynamical equations of massive LED. In any event, a 'state at an instant t' thus 
means the specification of the complete set of dynamical variables on the leaf with index t 
of a compatible foliation of spacetime. 

To find out what constitutes such a complete set of dynamical variables, we notice that 
our dynamical equations are of first order in the proper-time derivative of each of the kine- 
matical particle variables q, u, ilg, and of first order in the Lorentz time derivative of 
the fields F. A dynamical state in our massive LED will therefore be a set of the form 
{q(r), u{t), rtE^r), F( . )}, specified on a leaf St of a compatible foliation of spacetime in- 
dexed by t, satisfying various consistency conditions that constrain the values of q, u, Oe at 
T and the field F( . ) on the leaf. In particular, the four-velocity has to satisfy ||u|| = —1, 
and the anti-symmetric rank-two Euler tensor has to satisfy ^^^(t) ■ u{t) = 0. Furthermore, 
the Maxwell-Lorentz equations impose the constraint that F satisfy the three-divergence 
part of the Maxwell-Lorentz equations, which fixes the longitudinal part of F as a passive 
attribute of the charged particle. 

These are all local constraints imposed by relativity, which occur also in the traditional 
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limiting theories, i.e. the relativistic point charge mechanics in given external force fields, the 
mechanics of a point spin in given external torque fields, and the Maxwell electrodynamics 
with given source terms, in all of which cases the constraints are known to propagate in time 
if they are satisfied initially. This will be true also for our massive Lorentz electrodynamics. 

However, our covariant first-order dynamical equations are also mildly non-local, and 
this is going to be refiected in the compatibility conditions to be satisfied by the foliation 
of spacetime. Namely, since the first-order proper-time derivatives of u and J7e at proper- 
time r are given (implicitly) in terms of integrals of the fields over the particle's support 
in its instantaneous space slice of simultaneity, the transversal part (i.e. the dynamical 
degrees of freedom) of F(x) must be known for all x G supp (/e(||-||)) satisfying u{t) ■ (x — 
q(r)) = (recall that the passive longitudinal part of the field F(x) is not a dynamical 
degree of freedom). Hence, we say that a foliation of spacetime is compatible with a world- 
line r H^ q{t) of an extended Lorentz electron (i.e., a particle history) if each leaf of the 
foliation contains precisely one of the sets {u{t) ■ (x — q(r)) = 0} fl supp (/e(||-||)) (for 
simplicity, we assume that supp(/m) = supp(/c)). A compatible foliation sets the spacetime 
parameterization w.r.t. which the Maxwell-Lorentz equations must be solved in order to 
enable the evaluation of the Minkowski force and torque. 

A compatible foliation is of course not known explicitly until the particle world-line 
is known; it has to be determined along with the particle world-line and the fields by 
solving the second variation equations for the extrinsic curvature of the leaf S^ coupled 
with the constraint equations of Gauss-Codacci and Gauss in the complement of the sets 
X G supp (/e(ll-ll)) satisfying u{t) ■ (x — q(T)) = 0, which is a non-linear free boundary 
problem. Interestingly, these circumstances inject technical elements from general relativity 
into this special relativistic theory. 

Remark: No foliation is compatible with a particle history that is not ordered in the 
sense described earlier (see Fig.l and 2). 

The condition that the particle history be ordered imposes size restrictions on the accel- 
erating fields. Indeed, we easily find locally consistent data for which the particle's three- 
acceleration in its instantaneous Lorentz rest frame exceeds the value 1/R, where R is the 
radius of the particle. For such accelerations, part of the support of the particle in its instan- 
taneous space slice of simultaneity will be launched backward in Lorentz time rather than 
forward, and this motion does not generate an ordered history. This implies the following 
bound on the particle's initial three-acceleration in the Lorentz frame in which the particle 
is instantaneously at rest, 

\q\R<l, (8.1) 

which has to be translated into bounds on the accelerating field strengths in the instantaneous 
particle rest frame. We will soon see that these bounds on the field strengths are about 100 
times the particle's self-fields at it's surface which, happily, is way outside the physically 
reasonable range of data for a classical electrodynamics. Last not least, the condition that 
the gyrational speed for points on the particle's equator must be subluminal imposes a size 
restriction on i^E, viz. (recall that c = 1) 

||wE||i?<l. (8.2) 
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8.2 The Cauchy problem for the state 

We now formulate the initial value problem for the states. We begin by specifying the generic 
initial data on a convenient initial leaf. The compatible foliation itself is characterized 
subsequently by dynamical equations which together with the equations of massive LED 
(rewritten w.r.t. the compatible foliation) pose a nonlinear, implicit Cauchy problem for the 
joint evolution of the state and the leaf of the foliation on which the state is defined. 

8.2.1 Cauchy data and the initial leaf 

Since by a suitable spacetime translation we can always transform to a Lorentz frame in 
which the event of the particle's initial spacetime position is at the origin of that Lorentz 
frame, and since by a suitable Lorentz boost we can furthermore transform to another Lorentz 
frame with the same origin but in which the particle's initial three- velocity vanishes, and 
since finally we can choose the proper-time of the particle equal to zero at that initial event, 
the Cauchy data for the world-line r i— > q(r) can be chosen conveniently as 

q(0) = O, (8.3) 

u(0) = eo . (8.4) 

Notice that eo satisfies the constraint ||eo|| = —1 required of a four-velocity. Moreover, by 
at most a space rotation of this Lorentz frame we can achieve that vv^e(O), the dual vector 
to r2E(0), has only a single non- vanishing component, say the 3 component, We(0) = cUeo^s, 
with Weo satisfying (|8.2| ) but otherwise at our disposal. The Cauchy data for the particle's 
gyrograph r i— >■ ^IeIt) are then of the form 

riE(O) =cuEoei Ae2. (8.5) 

Clearly, f2E(0) is anti-symmetric, and since (ei A 62) ■ eo = 0, also of space-space type w.r.t. 
u(0). 

The Lorentz frame in which ( ^.3|) , (|8.4|) , and (|8.5|) holds will be referred to as the initial 
rest frame Ti^q. The initial leaf So of all of our foliations is now chosen to be the particle's 
afiine space slice of simultaneity in the initial rest frame JFlo, i.e. the set {eo • x = 0} on 
which the Lorentz time t (= a;°) = 0. 

Accordingly, the Cauchy data for F(x) are specified on So, 

F(x)|i=o = Fo(a;), (8.6) 

where Fo(a5) satisfies the three-divergence equations in the Maxwell-Lorentz equations re- 
stricted to Sq. To make this somewhat implicit characterization of Fq more explicit, we 
decompose Fq into its electric and magnetic components w.r.t. JFlq and group them together 
into a complex electromagnetic three- vector field on Sq, 



dcf 



Cfo=fo + ^fio, (8.7) 

whose real and imaginary part are, respectively, the electric (i.e. time-space) and magnetic 
(i.e. space-space) components of the field tensor Fo in Sq. The three-divergence equations 
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in Eq now read 

V-go{x)=4npoix), (8.8) 

where po is the time component of J(0), uniquely determined in terms of the initial values 
of the other state variables Qq, uq and Oeq- 

Remark: Through the presence of u in the right-hand side of ( |5.iJ| ) the (initial) constraint 
appears to be implicit. However, by direct computation it is readily verified that the u con- 
tribution to J drops out initially as a consequence of the sphericity of the charge distribution 
in the particle's instantaneous rest frame, which initially coincides with the Lorentz frame 

The initial fields are chosen in accordance with the asymptotic condition that Qo{x) -^ 
as \x\ — > cx), the real part as Sq{x) ~ —ex/\x\^ + o(|a;|~^), the imaginary part sufficiently 
fast for surface integrals at infinity involving Bq to vanish. 

8.2.2 Evolution equations for state and compatible foliation 



To each differentiable particle world-line starting with ( p.3|) and (|8.4|) and generating an 
ordered particle history (referred to by an index h), we have to associate a convenient, com- 
patible foliation of spacetime which contains Sq as initial leaf. Unfortunately, the standard 
foliation is generally not a compatible foliation, for notice that any world-line is the soul 
of the tubular set of events UT-^oSupp(/e(||x — q(r)||)), called the future world-tube of the 
particle, inside of which the foliation is necessarily given by the r-indexed affine segments 
{x|u(r) ■ (x — q(T)) = 0} which intersect with the standard leaf of Lorentz time t but are 
generally strict subsets of that leaf. However, one expects that it is possible to select a 
compatible foliation whose leafs, sufficiently far away from the world-line, will arbitrarily 
closely approximate the leafs of the standard foliation of the initial Lorentz frame, i.e. Eg 
parallelly transported in Lorentz time. By further choosing the t-parameterization of the 
compatible foliation to be Lorentz time at spatial infinity, by selecting the corresponding 
value of r through the identification g°(r) = t, and by parameterizing points in St by the 
points a? G So through following the integral curves of VgT/i, we obtain a foliation which 
nearly coincides with the standard foliation in most of spacetime, except for a small "wiggle" 
accompanying the world-line of the particle. Notice that the world-line is one such integral 
curve, so that in this parameterization the segment of the world-tube that belongs to Sj 
is always centered at the origin of Sq. Equipped with such a global parameterization of a 
compatible foliation, spacetime is diffeomorphic to the product manifold M x Sq, i.e. M x M'^, 
with spacetime metric 

ds^ = -ij\x,t)dt^+ Y, l''Kx,t)dx^dx^ (8.9) 

where ipix,t) = ||— VgT/i(x)|p is the lapse function, and the 7'^'(cc,t) are the components 
of the first fundamental form 7 of the foliation induced on S^ from the Minkowski metric. 

Such an extension of the family of affine segments of the future world-tube to the outside 
of the world-tube is found by solving the structure equations for a foliation of Minkowski 
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spacetime, supplemented with standard foliation asyniptotics at spatial infinity, and with 
appropriate boundary conditions along the perimeter of each segment in the tube which 
guarantee a differentiable transition at the segment's boundary. The boundary is determined 
by q(r) and u{t) (given /e, which as a model parameter does not vary during the dynamics). 
Since we have to determine the world-line along with the foliation outside the world-tube, 
our problem belongs to the category of the dynamical free-boundary problems. Notice that 
even if the world-line is analytic, the foliation is surely not analytic unless the particle is 
stationary, in which case one can choose the standard foliation of spacetime. 

By following [^, ^, we conclude that outside the world-tube the foliation is a solution 
of the system of evolution equations 

dtj = -2ipK (8.10) 

dtK = -V^ ® V^ip + ilj{Ric + KTtK-2K-K) (8.11) 

where K is the second fundamental form of the leaf of the foliation, i.e. the extrinsic curvature 
(actually defined by ( ^.101 ) if T^ were known), Ric is the Ricci tensor of S^, and V^ denotes 
covariant gradient w.r.t. 7. Equation ( |8.11|) is known as the second variation formula. In 
addition, 7, ip and K have to satisfy the constraints given by the Gauss-Codacci equation 

V^-K-V^TrK = (8.12) 

and the Gauss equation 

Ric-K-K + KTiK = (8.13) 

which propagate in time if satisfied on the initial leaf. This set of equations has to be 
supplemented by boundary conditions along the perimeter of the world-tube segment for ip 
and for K; and by asymptotic conditions as |a;| — > 00. We demand that the foliation be 
normal, i.e. that ip{x,t) — i> 1 as \x\ -^ 00, which means that at spatial infinity the foliation 
tends to the standard foliation of spacetime. Furthermore, a "gauge" condition is at our 
disposal, for instance the maximality condition 

TyK = 0. (8.14) 

The Maxwell-Lorentz equations for the field tensor F have to be rewritten w.r.t. the 
foliation in terms of the associated decomposition of F into electric and magnetic components 
(cf. [0), here conveniently grouped together as a complex electromagnetic three- vector field, 

g=£ + iB, (8.15) 

the initial decomposition being defined already earlier. The Maxwell-Lorentz evolution 
equations for Q then read p6i 



dtG{x,t) = -tV^x {^g){x,t)-A'K{tljj){x,t), (8.16) 
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where dt means first-order foliation time derivative (wliicli equals Lorentz time derivative far 
away from the particle), V-^x is the curl w.r.t. 7^^'^, and the electrical current density three- 
vector j is the Ej projection of J. Clearly, ( ^.16| ) is an explicit equation for the foliation 
time derivative of Q; however, note that j{x,t) is a (vector) functional not only of all the 
other state variables Og, q, u but in addition also of u. 

Beside the evolution equation ( ^.16|) , Q has to satisfy the three-divergence equations 



V^-g{x,t)=A7rp{x,t), (8.17) 

where p, the time component of J, is likewise a functional of all the other state variables 
Oe, q, u and in addition also of u. However, ( |8.17|) is already included in our set of conditions 
to be satisfied by a state, for notice that by taking the (three-) divergence of (|8.16|) it follows 
that a solution Q{x,t) of ( |8.16| ) for J(x) given in terms of consistent q(T) and f2E(r), 



automatically satisfies (|8.17| ) for allt > if it satisfies ( p.l7|) initially. 



The Maxwell-Lorentz equations are supplemented by the asymptotic condition that 
Gix) ^ as \x\ -^ 00, the real part as S{x) ~ —ex/\x\^ + o(|a3|~^), the imaginary part 
sufficiently fast for surface integrals at infinity involving B to vanish. 

Finally, following [^ one converts Q{x,t) into F on a given leaf, which inserted into the 
world-line and gyrograph equations ( [Tsl) and (|7.6| ) closes the system of evolution equations 
for the state in massive LED. 

In summary, the Cauchy problem for the state evolution turns out to bear a surprising 
resemblance to the general relativistic evolution problem of a black hole interacting with 



electromagnetic radiation as described in |^6[, and to the evolution problem of the vacuum 



Einstein equations [^. This gives the equations of massive LED an unexpected additional 
appeal as a somewhat simpler "toy problem" for these more formidable type of problems of 
general relativity. At the same time, the enormous progress recently made in the rigorous 



treatment of the Einstein equations ||35| gives hope that the complete control of the Cauchy 
problem for the state in massive LED can be achieved along similar lines of analysis. 

We shall have recourse to the notion of the state in massive LED in our section about 
scattering. In particular, the invariants stated earlier have to be converted into conservation 
laws for the evolution w.r.t. a compatible foliation. In the remainder of this section now we 
will discuss an alternate, global way of looking at massive LED, including an algorithm for 
the actual computation of solutions to any degree of precision. 

8.3 Nearly free global evolutions 

If, instead of evolving the state locally in time, we look at the putatively harder problem of 
global particle histories and field decorations of spacetime, one gains however the advantage 
that one can work with the standard foliation. The price to pay, of having to solve the global 
problem, will be somewhat lessened by noting that for all practical purposes of computation, 
one will need only an approximate solution. We restrict ourselves to nearly free global future 
histories and decorations, with data prescribed at t = r = on Eq, as before. 

Nearly free evolutions start from Cauchy data which are sufficiently close to a stationary 
state. For such evolutions, the system of fully implicit massive LED evolution equations can 
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be converted into an equivalent, regular quasi- explicit first-order system. We then explain 
that this quasi-explicit system of equations can be treated via a generalized Picard iteration.^ 

8.3.1 Maxwell— Lorentz equations revisited 



We now work with the standard foliation, so that (|8.15| ) is the decomposition of F into 
electric and magnetic components w.r.t. this standard foliation. The Maxwell-Lorentz 
evolution equations for Q now read simply as 

dtG{x, t) = -iV X g{x, t) - A7TJ{x, t) , (8.18) 

where dt means first-order Lorentz time derivative, Vx is the conventional curl, and the 
electrical current density three-vector j is the space component of J. Clearly, ( p.l8| ) is an 
explicit equation for the Lorentz time derivative of Q. Once again, j{x,t) is a (vector) 

functional not only of all the other state variables fi^, Q, u but in addition also of u. 
Beside the evolution equation ( ^.181) , Q has to satisfy the three-divergence equations 



V -Gix^t) =ATrp{x,t), (8.19) 

where p, the time component of J, is likewise a functional of all the other state variables 



fig, q, u and in addition also of u. As explained earlier, ( p.l9| ) will be satisfied if it is satisfied 
by a initial state. 

The Maxwell-Lorentz equations are supplemented by the asymptotic condition that 
Q{x) ^ as \x\ -^ cxD, the real part as S{x) ~ —ex/\x\^ + o(|a;|~^), the imaginary part 
sufficiently fast for surface integrals at infinity involving B to vanish. 

8.3.2 Gyrograph equations revisited 

The gyrograph equation ([7^ ), when viewed as an equation for ilg, does not have the desired 

o o 

format in the first-order proper-time derivative of 0-^ because Sb involves 0^ in a rather 
complicated manner. However, it is possible and also convenient to view (|7.6|) as an equation 
for Sb. To see that ( |7^ ) can indeed be viewed as an equation for Sb, we recall that we 
required /m to generate a strictly positive moment of inertia Xb. As a consequence, the map 
||we|| I— > Xb(||wE||) is strictly positive, increasing, and strictly convex for ||we|| G [0,1/R). 
This implies that we can invert Sb = Xb(|| WglDfiE to get the Euler angular velocity tensor 
Ob uniquely in terms of the bare spin tensor Sb, viz. fl^ = X^^ (||Sb||)Sb, where X^^ (^) = 
^^"'^(Xbid)~"^(^). With Oe understood in this way as unique (vector) functional of Sb, the 

^ The terminology of "explicit, quasi-explicit, and fully implicit first-order system of equations," is readily 
explained at hand of just two coupled dynamical variables S,{t),C,{t). (A) a pair of equations ^ = fi{^,0: 
( = /2(Cj C) is called an expHcit first-order system. (B) a pair £_ = gi{£,, C, C), C — .92(C, £,, C) which is not also 
of the form (A) is called a regular quasi-explicit first-order system. Thus, in a regular quasi-explicit system 
each dynamical variable satisfies an explicit equation conditioned on the histories of the other variables. (C) 
a pair -^1(^,^7 CjC) = Oi P2{£,,^,CX) = which is not also of the form (B) or (A) is called a fully implicit 
first-order system. The class of fully implicit systems contains the singular quasi-explicit systems, which are 
of the form hi{^, C, (■)C — 5i(Cj Ci C); ^2(^7 Ci C)C = 32(^1 d C) with coefficient functions hi and /i2 that vanish 
somewhere in the range of their arguments. 
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Minkowski torque t given in (7^) is now a (tensor) functional of Sb, F, q, u, u, so that (JT] 
as equation for Sb, 



-^Sb — — [i^FW, SbJ _ + t , 



dr 



^.20) 



is already in the desired format regarding Sb- Notice that not only t but also Op^ still 
depends on u. 

8.3.3 World-line equations revisited 

We now come to the world-line equation (|7.3|) , which is of second order in the proper- 
time derivatives of q. Since we can treat q and u as independent degrees of freedom, this 
second-order equation is equivalent to two first-order equations; namely, ( [7.3|) , now viewed 
as equation for u, together with 



dq 
d^ 



u 



(8.21) 



for q. While ( ^.211) already is in the desired format, (|7.3|) is not. In fact, u enters through p 



in the l.h.s. of (|7.3|) and through the tensor Op^ = u A u in the r.h.s. of (|7.3|) , so that (|7.3| ) 
is fully implicit first-order no matter whether viewed as equation for u or for p. Switching 
to p as dynamical variable therefore does not seem to offer an advantage in the conversion 
of ( [7.3| ) into the desired format. 

Continuing thus to handle (|7.3| ) as an equation for u, we first sort out all terms involving 
u. Beginning with the r.h.s. of the world-line equation ( [7.3| ), we notice that f can be 
decomposed into a sum of two terms, one of which is independent of u and the other one 
linear in u, 



f=[ F{q + x)-{u-n^- x)/e(||x||)5(u • x)d^x 

F(q + x).u(u-x)/e(||x||)5(u-x)d^x. 
i 

The l.h.s. of (|7.3| ) can be decomposed in a similar fashion. For p = M ■ u we have 



p = M u+ M u, 

with M given in ( [4.2|) , and with 



(8.22) 



(8.23) 



M = Mbg 



+ 



X ® X, [F(q + x), Op] J ^/e(||x||)<5(u ■ x)d^x 

X ® X, [u ■ VgF(q + x), flp] J ^/e(||x||)(5(u ■ x)d^x 

X ® X, [F(q + x), flp] 1 /e(||x||)(u ■ x)u ■ Vg5(u ■ x)d^x . 



(8.24) 



26 



The tensor M is independent of u, and so are the first three terms of M, while the last term 
in M is linear in u. In conclusion, the world-line equation can therefore be written as 



M u = f , 

where the four-force f depends only on q, u, J^e, ^e, F, and dff, thus 
f = -M^u + [ F(q + x) ■ (u - f^E ■ x)/e(||x||)5(u ■ x)d^x 

+ / x(x ■ Oe • (u ■ VgF(q + x)) • u)/e(||x||)5(u ■ x)d^x 

jRi.3 

+ f x{x- [F{q + x),n,] ■ u)U\\x\\)6iu- x)d'x, 

jRi.3 

and the pseudo-inertia tensor M only on q, u, fig, F, and dfP, thus 



(8.25) 



(8.26) 



M 



Mbg 



x®x, [F(q + x),riJ^ /e(||x||)(5(u-x)d^x 
/ x(g)x(x- r^E • F(q + x) • u)/c(||x|| )u- Vg5(u- x) d'^x 

JRi.3 

/ F(q + x) ■ u(g)x/c( ||x||)5(u-x)d'^x. 

JR1.3 



(8.27) 



Remark: The first three terms in (|8.27| ) are manifestly symmetric, but the fourth one is 
not and, as far as we can see, also cannot be symmetrized by adding a term whose action on 
u yields zero. Although this brings some technical inconveniences for concrete computations, 
it does not interfere with the remaining steps in our analysis. 

At this point, the system of evolution equations ( ^.18|) , (|8.20|) , (|8.21|) , and (|8.25|) is already 
a quasi-explicit first-order system; however the tensor M could be singular (as an operator 
acting by matrix multiplication (w.r.t. the Minkowski metric) on four- vectors) . Our next 
step is to show that in typical situations M is not singular. We will establish the invertibility 
of M for admissible initial data in a large neighborhood of those data that correspond to 
the stationary solutions with empirical electron characteristics that covers essentially all 
physically sensible situations. 

We begin with the observation that, since we have chosen a Lorentz frame for the initial 
value problem which coincides with the particle's initial rest frame, we can in particular 
choose initial data such that the particle and the fields remain stationary for all time. The 
stationary electric Coulomb field and the stationary gyro-magnetic field field, denoted here 
by an index st, are defined by the equations 



-iV X Gst{x) = 47r/c(|a;|) lJe x x , 
with cJe = <-fEO constant, and 

V-gst{x) = 47rmx\) , 



(8.28) 



(8.29) 
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together with the condition that Gst{^) — > as |a3| — i> oo. For a stationary situation in 
which the charge, gyro-magnetic moment, and renormahzed mass are matched to those of 
the empirical electron, with rr^b > but small (details are given in section 10.1), the first 
and third field integrals in (|8.27|) , evaluated with the corresponding Fst, are of magnitude a 
(Sommerfeld's fine structure constant) relative to the gyrational bare term, while the second 



field integral in ( 8.27 ), evaluated with the corresponding F^t, vanishes. The field integrals in 
( ^.27| ) are therefore only small perturbations (in the mathematical sense of operator theory) 
to the gyrational bare term A^bg in ( p.27| ), which in turn is invertible because it acts just 
as multiplication by the strictly positive coefficient A^b (~ f^e, the empirical electron mass) 
— recall that the metric tensor g acts as identity on four-vectors. Hence, we conclude that 
for a stationary, empirical data-honoring situation, M is invertible. 

Having established the invertibility of such a stationary M, we now turn to the discussion 
of initial data which do not launch a stationary solution. Since a < 10~^, in such a non- 
stationary situation the tensor M will remain invertible roughly up to the point where the 



field integrals in ( p.27| ) become comparable to the gyrational bare term. Inspection of the 



field integrals reveals that this limit is reached (a) for applied initial field strengths about 
100 times the surface field strength of the stationary electron, i.e. about 10®V/m, or (b) 
for values of — zV x Qoi^) ~ 47r/e(|a;|) uj^o x x about 10^ times the surface field strength of 
the stationary electron per radial size of the electron. Interestingly, field configurations that 
strong are precisely those that saturate the admissibility condition (|8.1|) . While therefore it 
is to be expected that a proof, if at all possible, of the general invertibility of M for all a 
priori admissible initial data will be a technically hard problem, from a physical perspective 
it is idle to worry about field strengths that strong within the Lorentz model. These critical 
field strengths are way outside the reasonable range of applied field strengths for which one 
can hope that the Lorentz model makes physically sensible predictions; whence, we will not 
address the issue of the critical field strengths any further and simply conclude with the 
observation that for essentially all physically acceptable initial data the tensor M is a small 
perturbation of A^bg and therefore invertible. 

We may now multiply ( p.25| ) by the inverse of M to obtain the first-order evolution 
equation 

^ = M'f, (8.30) 

dr 



valid surely initially and, if a continuously differentiable solution u of ( 8.25 ) exists, valid 
also for at least a short ensuing time interval during which M will continue to be invertible, 
then, literally by continuity. Conversely, if (^.30) has a continuously differentiable solution 



u, then as long as M remains invertible this solution will also solve (|8.25| ). Equation ( p.30|) 
is called the hodograph equation for u. It is of the desired format in u. 

8.3.4 Consistency 

Before we address the question of how to solve our quasi-explicit system of evolution equa- 
tions, a remark is in order regarding the following subtle point. Since u • u = 0, by taking 

-1 



the inner product of ( p.30| ) we see that we must have u ■ M ■ f = 0. Unfortunately, this 
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last equation is not an automatic identity (as is, for instance, u • F ■ u = 0, which follows just 
by the anti-symmetry of F, irrespective of whether u and F are dynamically related or not), 

and this may cause some concern. However, it is only necessary that u ■ M • f vanishes 
identically for solutions to the dynamical equations, and this can be shown to be the case. 
In fact, it holds for any solution of ( ^.301) with given fields F and gyrograph fig conditioned 
on r I— > OE(r) being a solution of ( [7.61) . In brief, the apparent difficulty traces back to 
the fact that Nodvik's Minkowski force f given in (|7.5|) is generally not four-orthogonal to 
u |T^ because it couples the spin-generated electromagnetic (dipole and higher multipole) 
moments to the inhomogeneities of the electromagnetic fields over the size of the particle. 
We have 

fu = - f {x-q)-n,- F(x)/e(||x -q\\)6(u-{x- q)) d'x-u, (8.31) 

which vanishes identically (irrespective of F) only for a so-called 'particle without spin' for 
which fig = (more on particle models without spin in Appendices A. 2 and A. 3). This non- 
orthogonality of four- velocity and Minkowski force is countered on the left side of ( |7.3| ) by the 
spin-orbit coupling term M^ in M, which creates a generally anisotropic translational inertia 
of the gyrating charged particle so that the Minkowski momentum p and the four-velocity 
u are not anymore parallel to each other but related by the genuine tensor proportionality 
([r.4|). By adapting the corresponding computations of Nodvik |1^ to the present model, it 



can now be shown that the equation obtained by taking the inner product of (|7.3| ) with u is 
identically satisfied conditioned on r i-^ ^e{t) being a solution of the gyrograph equations 

( [r.6|) , ( [7.7|) . The correspondingly conditioned identity u ■ M ■ f = now follows from the 
fact that ( |8.30| ) is equivalent to (|7.3|). 

8.3.5 Iterative treatment 

Having succeeded in writing the whole evolution problem as a regular quasi-explicit first- 
order system ( p.l8|) , ( |8.20| ), (|8.21|) , and ( p.30|) , in the final step in our local analysis of the 



initial value problem we now explain that this regular quasi-explicit first-order system can be 
treated by a direct generalization of the familiar Picard iteration scheme. This then verifies 
that a unique regular solution exists and, above all, shows how it can actually be computed. 
Thus, we consider the coupled sequences of approximate world-lines {q^'^\T)}'^^Q, r ^ 0; 
approximate hodographs {u'^"^(r)}^Q, r ^ 0; approximate gyrographs {CIe {t)}'^=q, t ^ 0; 
and approximate electromagnetic fields {Q^'^\x,t)}'^^Q, with a; G M^ and t ^ 0, jointly 
defined by 

— q("+i) = u(") , (8.32) 

dr 

Au("+i) =(M"Y"^f^"\ (8.33) 

dr ^ ^ 

l_S,(-+i) = _ [n^Jn)^ sj")] _ + t(") , (8.34) 
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and 

9i^("+i)(a;,t) = -iV^ X (V^C7)(")(a;,t) -47r(^j)(")(a;,t) , (8.35) 

together with the initial-time constraint 

V ■ ^("+') {x,0)= 47rp(") (a;, 0) , (8.36) 

starting the iteration with functions that are constant in time, 

(8.37) 
(8.38) 
(8.39) 
(8.40) 

and with correspondingly vanishing time derivatives of the starting data, 

54^(o)(x) = (8.41) 

q(o) = (8.42) 

uW = (8.43) 

Sb(°)=0. (8.44) 

This set of equations determines the first-order time derivatives of the maps t \-^ F^"^ ', 
T I— »• q^"-'^^\ T I— > u^"'~^^\ and r i-^ Sb , in terms of the order n maps t i— *■ ^'•"\ r \-^ q''"\ 
r I— > u^"'\ and r i— > Sb , and the first-order time derivatives of the order n maps. The 
maps t \-^ F'-"''' ' , t \-^ q^^~^^\ r \-^ u^'"'^^\ and r \-^ Sb themselves are now obtained 

by simply integrating the whole system with respect to the respective time variables. This 
is done for the Maxwell-Lorentz equations ( p.35| ), ( p.36| ) by the standard Laplace transform 



technique, thus picking up the consistent initial data Qq (i.e. Fq) posed at t = through 
the familiar integration by parts on dtG- Equations ( |8.32| ), ( ^.33| ), and ( |8.34| ) are simply 



integrated w.r.t. the proper-time from r = to r, thereby picking up the initial data q(0), 
u(0), and Sb(0) (viz. r2E(0)) from the integration of the explicit proper-time derivatives in 
the respective left-hand sides. This completes the n + 1^* iteration cycle. Starting with the 
prescribed order zero maps ( |8.4U| ) and their derivatives ( ^.44] ), one can thus systematically 



proceed to higher order iterates. The upshot is an algorithm for the practical computation 
of the solution, provided the iteration converges. 

Remark: It is important to notice that "practical computation" means to be able to 
obtain an approximate solution to any degree of precision in some finite region of spacetime. 
To reach this goal it is not necessary to actually compute the global iterands in spacetime 
(which would be impossible practically). A finite number of iterations on a nested, decreasing 
sequence of spacetime regions, starting with a sufficiently big initial region, will suffice. 
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8.3.6 Some remarks on rigorous results 

The convergence of our algorithm is equivalent to a contraction principle for the fixed-point 
map in the global future particle histories starting with the specified Cauchy data, their 
time derivatives (treated as independent functions in the iteration and in the fixed point 
map), and the field decorations of spacetime in the future of So, starting with the initial 
field configuration on Sq. Proving the contraction principle establishes the convergence of 
the iteration to a unique global solution in the future (and under similar conditions also the 
past) of the initial time for those admissible Cauchy data for which the tensor M remains a 
small perturbation of TWbg and therefore invertible. 

So far, for the model in which the bare rest mass and the charge are both concentrated 
on the same sphere, and provided the ratio of electrostatic to bare rest mass is smaller than 
1, we proved global existence and uniqueness for the iteration algorithm in the special case 
that only the particle's spin degrees of freedom participate in the particle + field dynamics. 
In the further specialization that the axis of rotation remains fixed, we were able to prove 
exponentially fast approach to a stationary gyration state for the field-particle system, con- 
vergence understood on families of nested compact sets, superposed on which is an outgoing 



field of electromagnetic radiation. For the proofs of these results, see |]33[ . 

For more general evolutions, convergence has not yet been established. The dynamical 
equations of our massive LED may spoil the global existence of the dynamics in several ways. 
First, the equatorial rotation speed has to remain subluminal for a particle with finite bare 
rest mass m-^. (As we will see later on, the situation is different in the limit m\^ — > where the 
rotation speed can reach the speed of light, or actually is locked in at that speed). Second, 
the acceleration of the particle has to remain below a critical threshold (mentioned above) to 
guarantee an ordered particle history. Third, but only relevant to the multi-particle systems 
not considered here, any two particles have to have a minimum distance between them to 
avoid overlapping cutoffs. We now comment on the first two problems. 

The first problem can be controlled by an appropriate a priori choice of the bare rest mass 
distribution. It suffices to discuss two representative choices, volume and surface inertia. 
They differ in their limiting behavior when ||vv^e|| -^ 1/R. The volume choice for f^ is 
representative for the larger class of continuous functions /m € C°([0,i?]), in all of which 
cases lim||vt^jj||^i//j A^bdl WeII) < cxd. In such a case it takes only a finite amount of gyrational 
energy to bring an equatorial point to the speed of light. This signals that the choice of a 
continuous /m will create a mathematically delicate problem when it comes to the question 
whether the gyrational dynamics remains subluminal, and so one may expect that this choice 
may easily lead to a singularity in finite time. In contrast, for the surface inertia choice for 
/in, the function A^b is not only increasing and convex, but satisfies A^bdl^Ell) — > cc as 
II We II — >■ 1/R. In this case it would require an infinite amount of gyrational energy to bring 
an equatorial point to the speed of light, so that it is to be expected that the dynamics will 
not lead to a singularity in finite time. These considerations show that surface inertia is 
preferable, although we do not claim that this requirement is necessary for global or even 
local existence. 

The second problem is much harder to control a priori. For a moving particle one now 
has to rule out that constructive interference of electromagnetic radiation fields leads to a 
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violation of the acceleration limit (in our Lorentz frame) 

l«l<7^. (8.45) 



where R is the radius of the particle and 7 = 1/v 1 — q^- At the initial instant t = r = 0, 
the particle is at rest and the estimate ( ^.451 ) reduces to (p.ll). 



Notice finally that the smallness condition that we need in our proof in [Q is surely 
going to be violated if one lets the bare rest mass flow to zero. The control of this situation 
apparently requires completely different techniques. 

Due to all these technical difficulties, the global existence and uniqueness of regular 
solutions, at least for the physically relevant parameter range of the model, still awaits its 
verification. In the following sections we proceed under the assumption that the dynamics 
of the field-particle system is globally well-defined. Put differently, whatever is concluded in 
those sections applies only to those situations for which a unique solution exists for all time. 

9 Scattering 

Having established that the dynamical state in our single-particle massive LED in spacetime 
M}'^ is given by the transversal part of Q (viz., F) on a spacelike hypersurface of spacetime 
containing the particle support, together with the particle variables q, u, and fig (or Sb) 
coincidental with the particle support in the spacelike hypersurface, we now turn the im- 
portant question which of these states are scattering states, and in which sense. Loosely 
speaking, a scattering state is a state whose backward evolution tends asymptotically in 
the infinite past to a superposition of free evolutions of propagating (in space) field and 
particle degrees of freedom, and similarly its forward evolution tends asymptotically in the 
infinite future to a(nother) superposition of free evolutions of propagating field and particle 
degrees of freedom. The scattering problem is concerned with establishing the existence of 
the scattering states, with their unique identification, and further with their classification 
according as to how their asymptotically free evolution in the infinite past is connected with 
their asymptotically free evolution in the infinite future. 

We will limit our discussion to situations in which the conserved quantities P, L, and s^ 
are finite. Moreover, we take a foliation in which P = (M, 0), i.e. a center of mass frame for 
the field-particle system. We assume the center of mass is at the origin of the initial leaf. 

We expect that in our massive LED there are basically three categories of states, namely 
the non-propagating (a.k.a. bound) states, the scattering states, and the 'catastrophic 
states'. For the catastrophic states the evolution either terminates at a finite time, or blows 
up in infinity time, or suffers some other pathological behavior like loss of uniqueness. As 
explained in the previous section, the global existence and uniqueness problem is largely 
unsettled so that we will not be able to rule out catastrophic states nor characterize them as 
physically irrelevant (although we expect the latter to be true). However, the special global 
existence and uniqueness result in [Q shows that for a certain class of initial conditions 
the evolution does exist globally, and these states are in fact scattering states. We shall 
show in this section that the scattering states have the remarkable, and in fact physically 
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indispensable character of soliton dynamics for the renormahzed particle degrees of freedom. 
As for the bound states, we readily characterize the stationary states, but we have no control 
yet on the question of whether periodic (or perhaps quasi-periodic) non-propagating states 
exist as well. 

This section on scattering is mainly included for the purpose of demonstrating the im- 
portant solitonic character of the renormahzed particle dynamics in scattering situations. 
As such, we will necessarily be brief on the other scattering issues. 



As a general technical reference to scattering theory we recommend [0, though most 
of the emphasis there is on quantum mechanical potential scattering; in this context, see 
also ^8|, ^ for subsequent technical breakthroughs, and |^ for a very clear exposition of 



the conceptual issues involved and further references. Various classical scattering results in 



electrodynamics are discussed in 26 



9.1 The bound states 

Massive LED describes the interacting dynamics of two a priori subsystems, the charged 
bare particle and the electromagnetic field, each of which is characterized by its own dy- 
namical degrees of freedom. With only a single charged bare particle interacting with the 
electromagnetic field in otherwise empty space, it could seem that there are no bound states. 
However, some of the particle's spin degrees of freedom couple strongly to some of the dy- 
namical degrees of freedom of the electromagnetic field, as a result of which bound states 
exist. The most elementary ones give rise to the notion of the renormalized particle. In 
addition, more subtle scenarios may possibly lead to dynamical bound states, on which we 
comment only briefly. 

9.1.1 The stationary bound states 

The most obvious bound states are stationary solutions with nonvanishing particle spin. 
Clearly we can choose the standard foliation of spacetime for these states, i.e. a Lorentz 
frame JFl in which the particle is at rest at the origin of the space slice Sq ~ S^; i.e. 
q(r) = for all r so that u{t) = Bq. The stationary states compatible with our conventions 
made in the Cauchy problem have been defined earlier; all other stationary bound states 
that satisfy our convention about the center of mass are now obtained by at most a rotation 
from these stationary states. 

For later convenience, we here list the solutions explicitly. The stationary particle in the 
center of mass frame is spinning with constant angular velocity tensor Og. The stationary 
charge- current density four-vector thus reads 

Jst{x) = (eo-f^K-x)/e(||x||). (9.1) 

We seek a stationary solution of the vector wave equation D A^j = Jst satisfying the Lorentz 
gauge condition and vanishing at spatial infinity. Recalling that — Oe • x = (0, We x cc), and 
introducing the time-space decompositions for the current density J = (p, jf) and electro- 
magnetic potential four- vector A = (0,-4), for the stationary Coulomb potential (j)st{x) we 
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find 



^st{X) 



;^/e(l^1)d 



'x' 



\X — X' 

and for the stationary vector potential A.sti'X:)-, 

Ast{x)=u}^^ / ^ fJ\x'\)A^x'. 

Jr3 \x -x'\ ^ 



(9.2) 



(9.3) 



The electric and magnetic fields are obtained in the usual manner as £st{x) = — V0st(a;) 
and Bst{x) = V X Ast{x), so that 



Qst{x,t) 



-'st 



X 



iV X Astix 



(9.4) 



Since we chose a spherical distribution /e(| . |) with compact support in R^, the Coulomb 
and the vector potential take a universal form outside the support of the particle, given by 



%t{x 



1 



\x\ 



Ast{x) = fix 



X 



\x\ 



ioT r > R, 



(9.5) 



where /x is the magnetic moment for the stationary lOj^, 



1 



t^ = ljl XX {uj^xx)f^{\x\)<rx. 



(9.6) 



Clearly, away from its support the stationary charged particle's electromagnetic signature is 
that of an electric point charge and a magnetic point dipole, as was to be anticipated. 

Finally, we remark that the static limit Iu^eI -^ for the stationary bound state is no 
longer a bound state of the field-particle system, in the sense that only the passive Coulomb 
field remains attached to the charged bare particle. 



9.1.2 The non-stationary bound states 

In various approximate, linearized versions of (massive) LED, time-periodic non-propagating 
field-particle solutions have been found. Time-periodic solutions do not exist if the Wiener 
condition is satisfied, viz. the Fourier transform of /e be strictly positive. However, notice 
that neither the surface density nor the volume density discussed in this paper satisfy the 
Wiener condition. In principle one should therefore be prepared for the possibility of periodic 
solutions in our massive LED. On the other hand, rigorous studies of simpler semi-relativistic 
models of a particle interacting with a scalar wave field show that time-periodic bound states 
of their linearized version are structurally unstable to the nonlinear perturbation of the linear 
dynamics ETl; see also |E2l. No such study has yet been carried out for our massive LED. 
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9.2 Scattering states 

In our center of mass initial leaf, we can divide the scattering states into two categories: 
(a) scattering of electromagnetic radiation off a non-moving particle, and (b) scattering of 
radiation and moving particle in an encounter. 

9.2.1 Scattering with a non-moving particle 

Scattering with a non-moving particle is considerably simpler because the family of space 
slices of simultaneity for the particle are just the standard foliation of spacetime, as for a 
stationary particle. The only dynamically active degree of freedom of the particle is its Euler 
tensor, or in the stationary Lorentz frame of our foliation, the angular velocity vector ujTi{t). 
The dynamical degrees of freedom of the electromagnetic field are restricted to be compatible 
with the condition that the particle sit still at the origin of space. For instance, large classes 
of axisymmetric field decorations of spacetime will be fine. This case has been studied in 
by us in rigor; details will be given in |Q. Here we summarize the main results, obtained 



under the condition that the ratio of the particle's electrostatic Coulomb energy to the bare 
rest mass is smaller than 1. 

The best understood situation is when the particle's axis of rotation does not change 
during the evolution, say cUe = ^e^s, so that Uj, is the only remaining dynamical particle 
degree of freedom. For all admissible initial conditions corresponding to these restrictions, we 
know that on families of nested compact sets the field-particle system converges exponentially 
fast to a stationary state which is a bound state if the conserved s^ 7^ initially and the static 
state if s^ = 0, superposed on which in either case is an outgoing field of electromagnetic 
radiation. The evolution of the electromagnetic field is thus given by the scattering formula 

G{x,t) '^±^ g,,{x) + e-^'^xGZi^), (9.7) 

where ^°^^ is a divergence-free electromagnetic field, orthogonal to all bound-state fields, 
uniquely determined by the initial state, characterizing the outgoing radiation. Similarly 

g*{x,t) '^^ G*Ax) + e^'^>^g':^^*{x) (9.8) 

where * denotes complex conjugate. We remark that Qst is the same stationary bound state 
in both formulas, which is a simple illustration of the soliton dynamics of the renormalized 
particle in this special situation. 

In the general non-moving particle case where the axis of oJe is not fixed during the 
evolution, we know that every admissible initial state compatible with the constraints evolves 
uniquely into the future (respectively, past), but we do not yet know that on families of nested 
compact sets the field-particle system converges (exponentially fast) to a stationary state. 
However, we do conjecture that this will be the case. Assuming this to be true we were able 
to show that for large times the evolution of the electromagnetic field satisfies the scattering 
formulas 

G{x,t) *^°° C?r(^) + e-'*^^«(a.), (9.9) 
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and 



g*ix,t) '^^ ^rN + e^*^^a*N, (9.10) 



where now ^^"^ and Q°l^ are generally not the same stationary bound state; however, they 
differ by at most a space rotation, and this once again illustrates the soliton dynamics of the 
renormalized particle. 

In both situations, the explicit characterization of the scattering operator from the "in" 
states to the "out" states has not yet been worked out. 

9.2.2 Scattering with a moving particle 

Scattering with a moving particle can mean any scattering process in which the particle 
moves for at least a finite interval of time. In that case we need to work with a nonstandard 
foliation, as explained in the section on the Cauchy problem for the state. Needless to say, so 
far this is the least explored case. The explicit characterization of the scattering states, and 
of the scattering operator from the "in" states to the "out" states still have to be worked 
out. 

However, although it is rigorously unknown which initial states are scattering states for 
which the particle moves, a priori speaking there are several types of scattering processes 
to be distinguished: (A) the particle moves only a finite distance, associated with which 
are bound states in the distant past and future — this is a perturbation of the previous 
"scattering of radiation off a bound state;" (B) the particle is at rest only in the past 
(respectively, only in the future) but moves freely in the asymptotic future (respectively, 
past); (C) the particle moves freely in the asymptotic past and future. Furthermore, one 
can give the following partial characterization of the "in" and "out" states. 

In the situations (A) and the parts of (B) associated with a stationary bound state, the 
scattering asymptotics differs from the one discussed in the previous subsection by at most 
a translation of the stationary bound state. 

In the situations (C) and the parts of (B) not associated with a stationary bound state, 
the scattering asymptotics in the specified foliation is given as follows. As for the asymp- 
totic particle motion, all the "in" states {q™, u™, fiE™}, as well as all the "out" states 
|qOut^ jjout^ Oe°"*}, are obtained from the stationary particle states considered above by the 
action of an element of the Poincare group, more explicitly a combination of a translation, 
a rotation, and a Lorentz boost. The free evolution of a particle is obtained by applying the 
Lorentz-time propagation operator to these states, or expressed in terms of proper-time r 
(recall that q^{T) = t), 

q(r) "=±?° q°"* + u°"V (9.11) 

u{t) "^^ u°"* (9.12) 

OE(r) "=i°° Oe""* (9.13) 
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and correspondingly for r — i> — oo using the "in" states. Associated with these free asymp- 
totic evolution of the particle are the co-moving electromagnetic fields. Sufficiently far away 
from the particle these are just co-moving in the standard foliation to which the compatible 
foliation is asymptotically close except in a neighborhood of the world-tube. Their elec- 
tromagnetic "in" and "out" states are parameterized by the space part x°^^ of q°"*, the 
constant velocity v°"* in u"*^* = (7,7v)°"*, and the angular velocity vector oJe""* correspond- 
ing to JIe""*- To have a handy abbreviation for this parameter set, we choose the index com 
for 'co-moving.' We set 



^) ^^ ; + (1 - v') '^ ^ ^"^^ , (9.14) 






with 



v = v/\v\, y, = ix-x^^')-v, y^ = x-x°^'-y,v, (9.16) 



and where v = t;°"* and /x = /x°"*. After these preparations, we can give the asymptotic 
future evolution of the electromagnetic field for such a scattering state as 

Gi^,t) '^^ -v(e--°"^-^0°tN)-(5t-^Vx)(e---'v^°l(a.)) 

+ e-*^x«N (9.17) 

Similar formulas hold for the asymptotic past evolution in terms of the "in" potentials. 

9.3 Soliton dynamics 

While the rigorous treatment of the scattering problem of massive LED is still in its infancy, 
in all instances of scattering a remarkable conclusion about the gyrational particle degrees 
of freedom can be made. Naively one would expect that generally the orientation and the 
norm of the angular velocity of the particle would change during the scattering process — 
and with it the rotational and renormalized mass of the particle, its spin and magnetic 
moment. However, precisely this does not happen. Any scattering process connects two 
freely evolving particle states with identical renormalized mass and identical magnitudes of 
spin and magnetic moment. The relativistic Lorentz electron, equipped with no n- vanishing 
positive bare inertia, thus has the dynamical characteristics of an elementary particle in the 
best sense one could have hoped for in a classical theory. 

By adapting the terminology of Spohn |^3| and collaborators, who considered charged 



particles without spin, we call the renormalized particle a spinning charge soliton. As re- 
markable as this soliton dynamics itself, is the fact that it is a simple joint consequence of the 
conservation of s^ = —(1/2) TrS ■ S together with the fact that S = S for a stationary 
bound state, and the invertibility of the map fig i— > S for such a state. Indeed, we need only 
consider a stationary bound state, for we know that the asymptotic future (past) evolution 
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of the scattered particle dressed with its co-moving fields (the renormalized particle) is just 
a translated and boosted such state. 

Explicitly, for a stationary bound state we have, in the Lorentz frame and in dual vector 
notation, 



1 — |cl»e X a;P 



frai\x\)d'x+ [ XxAstix)U\x\)d^X, (9.18) 

7k3 



the first term being the bare spin, the second the electromagnetic spin of the particle. For the 
bare spin s^, we already saw, in our section on the Cauchy problem, that it is parallel to uJe 
and the map Sb ^~^ ^e invertible. Using our formula (|9.3|) for the stationary vector potential 
of a bound state, we see by explicit computation that also the total spin s is parallel to cje, 
and the map s i^ We is invertible. We have s = I{\u:^\)u}^, with / > and map x i— * xl{x) 
increasing, vanishing for x = 0. Now, \s\ = \u}j,\I{\u}j^\), so that from the invertibility of 
X I— > xl{x) we conclude that |s| determines IcJeI uniquely in a stationary bound state. Thus, 
from the constancy of s^ we conclude that |a;E| in the asymptotic future state has the same 
value as in the asymptotic past state. The solitonic character of the renormalized particle is 
proved. 

10 The limit of vanishing bare rest mass 

We carry out the renormalization group analysis for a simple example where the various 
dynamical functions can be computed explicitly. Specifically, we here consider 'densities' /m 
and /e which are Dirac measures concentrated on the surface of a two-sphere of radius R, 
viz. 

;^/m(|x|) = —JM^l) = ^Si\x\ - R) . (10.1) 

The elementary charge e will not be tempered with and therefore not displayed as variable in 
the arguments of the functions below. Beside having set the charge — e of the model particle 
equal to that of the real electron, we also demand that its total model mass m matches the 
empirical electron mass, m = rrie, and that its magnetic moment |/x| matches the empirical 
electron magnetic moment, |/i| = /ig. 

Remark: In this and the next section we resort to units in which the speed of light c 
is displayed explicitly, to facilitate the comparison with the formulas used in physical data 
books. 

The renormalized particle spin will be derived in our model in the limit of vanishing bare 
rest mass, whence is not to be matched to a 'corresponding empirical value,' but rather 
compared. Naively speaking it is to be expected that in our classical theory with ultraviolet 
cutoff the spin-to-magnetic-moment relation 

with Lande factor g ^ gci = 1 should hold. Interestingly, we will find g ~ 2/3. 
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We remark that, without further elaboration, there remains a certain ambiguity in match- 
ing our stationary model data to empirical electron data inferred with the help of quantum 
mechanics from data of dynamical experiments, viz. mg in scattering experiments, and /ie 
in resonance frequency experiments in a Penning trap. The ultimate justification for these 
choices has to wait for some future dynamical analysis of our model. 

10.1 The limit m-^ -^ 0+ for the stationary electron 

By performing at most a space translation, we can assume that the stationary particle, and 
its self-fields, are centered the at the origin of the space slice M^ of the Lorentz frame in 
which the system is stationary. Since the time- component of the spin four-vector and the 
angular velocity four-vector vanish in our stationary frame, after at most a space rotation we 
can furthermore assume that the spin four- vector of the stationary bare particle, Sb, and of 
its self-fields, Sf, points along the 63 direction. Since 63 = (0, 63), we can canonically identify 
these spin four- vectors with the familiar bare and field spin three- vectors Sb? -Sf, respectively 
the angular velocity four- vector with the angular velocity three- vector lJe, which moreover 
can be written as Sb = s^oB^, S{ = Sfe3, and lJe = uj^e^, with cJe > 0. With these conventions, 
and with (|10.1D, the relevant formulae for the renormalization analysis are the following. 



Beside the electric monopole moment (i.e. the charge — e), the magnetic dipole moment 
of the empirical electron is the second electromagnetic characteristic measurable in a classical 
manner. The magnetic moment three- vector /x given in (|9.6|) (and re-installing a factor 1/c) 
evaluates to 



t^ 



/ {\x\Hd-x(^x)S{\x\-R)d^x -u;^ (10.3) 

-l-R^oo^es. (10.4) 



2c47ri?2 
le 
"3c 



Since it points in the 63 direction as well, we will write /x = fie^, with /i read off from (|10.4| ). 
The mass m of the renormalized particle is defined by 

mc^ = Wi, + W{ . (10.5) 

Here, Wh is the gyrational energy of the bare particle, defined earlier in ( [4.2|) and evaluated 
for general spherical /m in ( [4.3|) . Re-installing a factor c^ and inserting now (|10.1|) for /m in 
([4.3|), we obtain 

47r-K^ JR3 y/l - \u:^ X x\^/c^ 

1 ^ f ^eR\ 

= T^hC — — artanh . (10.7) 

UJeR \ c J 

Furthermore, Wf is the electromagnetic field energy of the gyrating charged particle, 

^i=^ [ (\^stix)\' + \B,tix)\')d'x (10.8) 



57r Jm.3 

1 e V 2 u^^R^ 

IH 

2R\ 9 c 



l + ^^^l- (10.9) 
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Our system of equations for the stationary renormalization flow is given by 

rrie = m{R, LJjs, m\^) , (10.10) 

(Xe=WR,CO^), (10.11) 



with the function 'm{R,uJi^^m\,) given through ( |10.5D , ( |10.7D , ( |10.9| ), and with the function 
|/x|(i?, cUe) given through ( |10.4] ). Explicitly, we have 



c , fuj^R\ leM / 2uj^^R^\ ,__, 

^.^^.-^artanh(^j+---(l + --^j ' ^^'-^'^ 

lx, = \-R^uj^. (10.13) 

6 C 

These two coupled equations determine a curve in (i?, ciJE,''^b)-space for which we will study 
the limit m-b — > 0. 

Equation ( |10.13| ) can immediately be solved for uj^ in terms of the radius R, 

o/^eC 1 
UJv. = O- 



e R^' 
Inserted into ( |10.12| ), this gives 

1 e / Ue 1 \ 1 e^ 1 / ul l\ 

"= = ».iS3-arta„h (3^-) + --,- (l + 2§^j , (10.14) 

which rewrites into a formula for the bare rest mass as function of R, 

rn, = rn,3!^i^^!hliyp^. (10.15) 

artanh ( 3— — 
V eR^ 

We are interested in the inverse function R{m\y). In particular, we want to know whether 
limmb^o R can be taken without encountering superluminal equatorial rotation speeds. 

We readily read off from (|10.15|) that the right side is well-defined only on the R- 
interval (3/ie/e, oo). From ( |10.13| ) and ( |10.19| ) it furthermore follows that IcUeI-R < c for 



R G {3iJ,e/e, oo), viz., no superluminal rotation speeds occur for R in the domain of defini- 
tion for the r.h.s. of ( |10.15| ). Furthermore, m\y{R) is a monotonically increasing function, so 
that the limit R \ 3/ic/e for m\^ will answer the question whether one can take the limit 
rrib — > for R. A simple inspection of ( |10.15 ) reveals that mb — >■ frtf, as i? — >■ oo and that 



mb — > as i? -^ 3/ic/e. Hence, by monotonicity there is a unique limit of vanishing bare 
rest mass for the radius of the renormalized particle, given by 

lim R = 3—. (10.16) 

mt^O e 

It also follows from our discussion that the equatorial rotation speed in this limit is precisely 
equal to the speed of light c. 
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The renorinalization flow analysis reveals that the Lorentz program of a "purely elec- 
tromagnetic electron" can be completed successfully, in the stationary setting at least. The 
resolution of the puzzlement is that beside the static electromagnetic field energy contribu- 
tions considered by Lorentz, the mass of the properly renormalized purely electromagnetic 
particle consists of an additional contribution which we call photonic, for it is the mass of 
an object without rest mass which rotates with equatorial speed equal to the speed of light. 
This contribution was missed by Lorentz who calculated without bare rest mass from the 
outset. 

It is instructive now to make use of the empirical relation for the electron magnetic 
moment, fig = {1 + a)fiji, where 

1 he , , 

/iB = ^ 10.17 

is the Bohr magneton and a = 0.001159652... the anomaly factor (see, e.g., |Q). Up to 
corrections of 0{a), equation (|10.16|) then reveals that the radius of the renormalized purely 
electromagnetic electron is 

^- = IRc, (10.18) 

e 2 

where 

i?e = -^ (10.19) 

rrieC 

is the electron's Compton length. This result is quite remarkable in many respects. It shows 
that our Lorentz-covariant LED in fact contains the correct ultraviolet cutoff length where 
heuristic discussions traditionally place the limit of applicability of a classical field theory, 
viz. at the order of the electron's Compton length Re, rather than at the (lamentably 
so-called) classical electron radius, given by 

^ci = ^ , (10.20) 

which gives the size of a 'purely electrostatic classical electron.' Recalling Sommerf eld's fine 
structure constant, given by 

e^ 1 

a = — ^ , (10.21) 

he 137.036' ^ ' 

we have the relation 

Ra = aRc (10.22) 

so that the true Lorentz electron is actually about (3/2a) ~ 200 times larger in diameter 
than predicted by a purely electrostatic calculation. 

We conclude this subsection with a discussion of the total spin, given by 

s = Sk + S{. (10.23) 
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Here, Sb is the spin three-vector of the bare particle, which for m\y > evaluates to 

s, = mb^ / —^^^^^^=6{\x\ - R) d'x (10.24) 

= ™.cfl( L^arta„h(^) -^^] e3, (10.25) 



which, by inserting ( |10.13| ) into (|10.25|) , with ( |10.17|) for fie, becomes 



-'^::il(^^^T^-- (11)41)- 

Using now ( |10.15|) , we get 

^ 3 l + i2R/3R,r ( _lRc.( IK\\ _ ^l^a^ .,0 27) 

The second spin term, Sf, is the stationary field spin three- vector associated with the rotating 
charged particle, which evaluates to 

^f = --T^ / ^ >< Astix)6i\x\ - R) d'x (10.28) 

— / XX (eJx)xBst(x))Sx (10.29) 



47rc 
2f 
9 c 



-—u^Res. (10.30) 



We write s = (sb + Sf)e3, with Sb and Sf given in ( |10.27|) and ( |10.30|) . 



We consider the spin coefficients Sb and Sf as function of R. Letting R -^ 1.5i?c, i-e. 
taking mb — ^ 0, for the bare spin we get, 

]ira sy. = h(l-l^a) , (10.31) 

while for the field spin, recalling that u)^R — > c as rrib — > 0, we immediately obtain from 
( |1030| ) that 

2 6^ 2 

lim Sf = = -ah. (10.32) 

mb^o * 9 c 9 ^ ^ 

Finally, adding Sf and s\, and taking the limit equals the sum of the individual limits, whence 
the renormalized limit spin magnitude is given by 

s,,„ = lim s = -/I ( 1 - —a ) . (10.33) 

nib^O 2 \ 21 I 



Put differently, from ( |10.33| ) we read off that the spin magnitude s^^ of the renormalized 
electron satisfies 

s_ ^ hi (10.34) 
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to high precision, with relative corrections of order a, compared to our naive expectation 
s ^ h, and compared to s ~ h/2 in a quantum mechanical framework. 

Curiously, furthermore, a closer look at the corrections reveals that in our LED an ex- 
pansion of the f^-factor in the Lande relation (|10.2| ) in terms of powers of a holds. We set 
go = 2/3, then 

^LED = 9o(l + ^a+ (^) a' + O {a^) j (10.35) 

which is surprisingly similar to the familiar expansion in QED |^, |ll 
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10.2 The limit rrih -^ 0+ for the dynamical electron 

We finally come the dynamical renormalization problem. In contrast to the stationary prob- 
lem, no rigorous results concerning the existence of this limit are available yet. However, 
it is of interest to inquire into the features expected for a renormalized dynamical limit of 
vanishing bare rest mass of massive LED. 

First of all, the matching of empirical electron data in the dynamical renormalization 
problem with scattering states has to be carried out in the asymptotic past (or future), for 
we want to "send in" the particle with the correct physical characteristic particle data. The 
incoming-data matching in the asymptotic past will fix R and ?7ib for the rest of the motion, 
since as parameters of the model these do not vary during the motion. 

This matching in the asymptotic past will also fix the value of I^eI in the asymptotic 
future, which is due to the soliton dynamics of the renormalized particle. During the motion, 
of course, \u:^\ will in general vary. However, our results on the scattering with a non-moving 
particle indicate that IoJeI approaches its asymptotic value extremely fast. Extrapolating 
from this observation, we expect that the soliton dynamics of the particle in massive LED 
in fact will try to keep the angular frequency \lo^\ of the particle as close as possible to the 
asymptotic value in the free motions. By the same token the gyrational mass (hence, also 
the renormalized mass, which is a small perturbation of the former in the stationary state) 
of the particle should remain (in the instantaneous rest frame) very close to the empirical 
mass during the motion. 

An important conclusion can then be drawn if one can show that the gyrational bare 
mass of the particle with rrih > stays uniformly close to the empirical mass as rrih —>■ 
for an rrib-indexed sequence of motions. Namely, in that case the equatorial gyration speed 
of the renormalized purely electromagnetic particle must be locked at the speed of light 
during the whole motion! For assuming to the contrary that the gyrational speed would 
remain bounded away from the speed of light for at least a finite fraction of the evolution 
(of course still close to it in the uniform sense stipulated) , it would then follow that the bare 
gyrational mass converges to zero, while the electromagnetic mass would contribute at best 
an amount of order a of the empirical mass. This is in contradiction to its hypothesized 
uniform proximity to the empirical mass. 
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We infer that in the hmit rri]^ — ^ the bare gyrational mass A^b is independent of o^g, 
hence a constant of the motion identical to what we have called the photonic mass of the 
stationary particle; explicitly, 

rripi, = mc ( 1 - —a j . (10.37) 

By the same reasoning, in the limit of vanishing bare rest mass the bare spin vector becomes 
the photonic spin, a four-vector with only two independent variables. It's dual tensor is given 
by 

Sph = ^n(^l-^«^S, (10.38) 

where S is anti-symmetric, of space-space type, has norm 1, satisfying S-u = 0. The inertial 
gyro-frequency fig and the spin tensor S are related by a simple proportionality, 

^ rripC^ „ 

Oe = ^S, 10.39 

n 

emphasizing the characteristic Euler angular frequency meC^/h explicitly. Similarly, fig = 
/tSph, with 

h^ 1 — ^a 

After these preparations, we can write down the putative dynamical equations of renor- 
malized LED. The world- line equation of renormalized LED reads 

Ap= f F(q + x)-(u-(f^™ + «:Sph)-x)/e(||x||)5(u.x)d4x, (10.41) 

where 

p{r) = Mre„(r) ■ u{r) (10.42) 

is the energy-momentum four-vector of the renormalized purely electromagnetic electron, 
with Mren giveu by 



M, 



rriphg- [ [x®x,[F(q + x),/s:Sph]J^/e(||x||)5(u-x)d4x (10.43) 



now the renormalization limit of the tensor mass M. The gyrograph equation of renormalized 
LED reads 

/ xA (g+u® u) ■F(q + x)- (u-(rJFw + /«Sph)-x)/e(||x||)5(u-x)d^x. 

(10.44) 
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The Maxwell-Lorentz equations for renorinalized LED have their famihar structure, but the 
charge- current density is now given by 

J(x)= f (u-(fi,w + /s:Sph)-(x-q)) /e(||x-q||)5(u-(x-q))dr. (10.45) 

The state in renormahzed LED differs from that of massive LED in that only the ori- 
entation but not the norm of Sph is a dynamical degree of freedom. This imposes some 



consistency conditions also on the initial conditions. We remark that ( |10.44|) does not au- 
tomatically preserve ||Sph|| for an arbitrary electromagnetic tensor field F and world-line 
T \—>- q. However, it is only mandatory that the norm of Sph be conserved for consistent joint 
solutions of the world-line, gyrograph, and field equations; the situation is similar to the one 
related to the conservation of ||u|| by the world-line equation. 

11 Summary and Outlook 

In the previous sections we have constructed the first viable relativistic Lorentz electrody- 
namics with positive bare rest mass (massive LED). We studied its renormalization flow limit 
to vanishing bare rest mass (renormahzed LED). This modern completion of Lorentz' pro- 
gram of a purely electromagnetic theory avoids the superluminal pitfalls of Lorentz' purely 
electromagnetic calculations. In the limit we obtain a particle with correct empirical mass, 
charge and magnetic moment, though with spin 3h/2. Although the bare rest mass van- 
ishes, the particle retains a "photonic mass" resulting from luminal equatorial motions, 
which actually accounts for most of the particle's mass while the traditional electromagnetic 
self-energies contribute merely a correction of order a. The size of the limit particle is of the 
order of the electron's Compton length, about 200 classical electron radii. 

The finite size of this ultraviolet cutoff in the 'purely electromagnetic' limit endows our 
LED with a mildly a-causal appearance. However, the ultraviolet cutoff is itself an individual 
but non-local event which cannot be subdivided into a collection of individual point events. 
The cutoff moves as a non-local unit in a Lorentz-covariant way, in fact accelerated in a 
causal manner — in the sense that fields outside the cutoff event do not influence its motion 
— by the field integrated over the cutoff at any given instant of its proper-time. It is perhaps 
surprising that such an instantaneous, non-local action is compatible with Lorentz covariance 
and causality, yet it is. 

The existence of this physically viable properly renormahzed LED reinforces some hopes 
that LED may ultimately have some bearing on the quest for a consistent QED, with LED 
or some improved version thereof (LED+ say) as its (semi-) classical limit. Recall that the 
perturbative series for QED [^ is most likely divergent [^] and hence does not give a 
mathematically acceptable definition of QED. Many experts nowadays are inclined to be- 
lieve that the constructive quantum field theory approach to QED [^] will not lead to a 
mathematically consistent theory either [^, though the failure of this approach to QED 
has not yet been conclusively demonstrated [Q. There is also a strong anticipation that a 



physically correct and mathematically well-defined relativistic quantum theory must involve 
also the strong and weak nuclear interactions, as well as gravity, and all this supersymmet- 
rically |^9|, ^ ^. It should however not be necessary to unify physics at the Planck scale 
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in order to obtain a consistent working theory for electromagnetic phenomena above the 
electron's Compton scale. Inspired by Kramer's program of 'non-relativistic QED' [^, ^ , 
which aims at the quantization of LED's semi-relativistic predecessor (see Appendix A. 3), 
and which is expected to describe the low energy limit of QED [^ , one may want to inves- 
tigate the possibility of deforming our LED into a 'quantized LED', perhaps along the ideas 
of Kontsevich [Q, and it would be interesting to see whether that would teach us something 
about QED proper. 

Such fundamental hopes aside, our LED is interesting in itself as a nonlinear relativistic 
microscopic theory which couples particle and field degrees of freedom in a self-consistent 
manner. While the rigorous control of its semi-relativistic predecessor has made significant 
advances in recent years due to the efforts of Spohn and collaborators and others (discussed 
in the Appendix), the extension of this rigorous control to our relativistic LED has just 
started ||3^. The most surprising new dynamical challenge is the need to construct a com- 
patible foliation of spacetime along the evolution of the state in LED by also solving the 
vacuum Einstein equations as a free-boundary problem coupled to the evolution of the state. 
Unexpectedly, massive and renormalized LED may therefore serve as playground for the 
harder problem of black hole dynamics in general relativity. 

Among the future projects, one is certainly the rigorous extraction of effective equations 
of motion for the particle in the adiabatic regime, where to lowest order one expects equations 
of the type discussed in pSJ pGJ , with radiation reaction-corrections showing in the next order 



of approximation. 

Another interesting project is the many particle problem, and in particular the rigorous 



microscopic foundation of relativistic Vlasov-Maxwell theory ||5^, ^ |5^, ^, which so far 
could not even be attempted for lack of a consistent microscopic theory of electromagnetic 
particles and fields. 

A Appendices 

Appendix A. 1 contains the more technical derivations of the Euler-Lagrange equations for 
our massive LED, Appendix A. 2 discusses the singular point mass limit in which our massive 
LED reduces to Nodvik's massive LED, and Appendix A. 3 provides supplementary material 
about the semi-relativistic LED. 

A.l The Euler-Lagrange equations of massive LED 

In this subsection, we first collect various formulas needed to carry out the variations of the 
action functional. Next we apply these formulas to derive the Euler-Lagrange equations for 
the action principle. 

A. 1.1 Independent variation of the world-line 

To a given world-line r i-^ q(r) we may add a small perturbation r i-^ Sq{T) to get a deformed 
particle world-line q(r) = (q + (5q)(r). Here, r is the proper-time of the unperturbed world- 
line. The proper-time along the perturbed world-line at the point q(r) will in general not 
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coincide with the value r, i.e. the proper-time of the unperturbed world-hne. To first order, 
the difference in the proper-times along the deformed and original world-line reads 

6{dT) = -u-6qdT. (A.l) 

o 

Accordingly, the derivative of 6q with respect to r, viz. 6q, is generally not the first-order 
variation in the four-velocity, 6u, which is defined by taking the derivative of q w.r.t. its 
proper-time, then computing the difference u — u to first order. As a function of r, the 

o 

first-order variation in the four-velocity is given by the space projection (w.r.t. u) of 6q, 

(5u= (g+u® u) -Sq, (A.2) 



where g is the metric tensor of Minkowski space, defined in ( ^^.121 ). Taking the inner product 



of (|A.2| ) with u, and recalling that g acts as identity on four-vectors, we find that u- Su = 0, 
as required by (u + 6u) ■ (u + 6u) = —1 in first order. The variation of fl-pw under an 
independent change 6q, the gyrograph being kept fixed, is given by |jl5[ : 



Srtpw = — u A Sq + ilpw u ■ Sq. (A. 3) 
Defining quasi-coordinates 

S'd = -uA5u, (A.4) 
equation (|A.3|) can be rewritten into 



d 
^fipw = -^S-d + [flpw, 5^]- + f^Fw" • Sq , (A.5) 

dr 

which proves convenient later on. 

Moreover, although the gyrograph is fixed during an independent variation of the world- 
line, the Fermi-Walker frame changes, and with it the tensor of the Euler gyration, for 
the gyrograph is defined w.r.t. the Fermi-Walker frame. The variation in Oe under an 
independent change 6q is given by [jl5| 



SQe. = -uA {flE -Su) + nEU -Sq, (A. 6) 

whence 

o o 

Sn = -uA {n -du) - uA 6q + nu -dq. (A. 7) 

A. 1.2 Independent variation of the gyrograph 

The orientation of the set of three spacelike vectors {§1, 62, §3} relative to the spacelike 
vectors {ei, 62, 63} can be specified by Euler angles 9, 0, ip. Independent changes 69, 6(j), 
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Sip of the particle's history of Euler angles, i.e. keeping its world-line fixed, give rise to the 
variation lITSll 



5nE = 56» + [rj,56»]_, (A.8) 

where the anti-symmetric tensor 50 is given by 

50= Y^ 5^^^e^®e^, (A.9) 

with {e^} the tetrad of the frame T^^ and 59^^ given by 



5e^^ = 5(j) + cosO 5ip 

502^ = cos (j) 59 + sin 9 sin 5ip 

59^^ = sin (j) 59 — sin 9 cos (p 5ip 

59o' =0 for z = 1, 2, 3 , 

cf. [0 for details. 

A. 1.3 Euler-Lagrange equations for the world-line 

We first turn to the action of the bare particle, which we rewrite as 



(A.IO) 



Aare= / L^i^^, u) dr , (A.ll) 

where 

Lb(0„ u) '^ - / V^l-||f^E-xf U{ \\x\\ ) 5{u ■ x) d^x (A.12) 

is the Lagrangian for the bare particle. We first define the volume S bounded by two space 
slices defined by u{t) ■ (x — q(r)) = at r = ri and r = T2. Its translation by q is denoted 
by Sq = {x — q^r) ; x G S}. Under an independent variation q 1-^ q + 5q the resulting 
first variation of ^bare reads 

5Aare = f (^ Tr (^ ■ 5n^ + ^J± . 5u - L^u ■ 5"q\ dr , (A.13) 

where the last term under the integral results from the change in proper time along the 
perturbed world-line, see equation (|A.1| ). As for the variation of Lb w.r.t. u, we get 



^^" - / y^l-llfl.-xf /m(||x||)9„5(u-x)d^x (A.14) 

•J =.a 



5u 



u f ^l-llf^E-xf /„,( ||x|| ) 5{u ■ x) d^x, (A.15) 
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where we used 

du6{u-x) = -x{u-Vg)6{u-x), (A.16) 

then integrated by parts, carried out the differentiations, and noticed that terms proportional 
to u ■ X in the integrand vanish upon integration against 5{u ■ x). But with the help of the 
four-orthogonahty of u and 6u we see that 

^■5u = 0. (A.17) 

ou 

As for the variation of Lb w.r.t. fig, we find that 



where Sb is the spin tensor of the bare particle, defined earlier in ( |6.12| ). Thus, we have 



(A.18) 



<5Aare = T Q Tr (Sb " 5^^) - Uu ■ 5q] dr . (A.19) 



Inserting now ( [A.6| ) for 50e, with ( [A.2| ) in place for (5u, then using Q,^ ■ u = once. 



next integrating by parts terms containing (5q, recalling that 5q = at the boundaries of 
integration, then using again Q,^u = and also noticing once again that terms proportional 
to u ■ X in the integrand vanish upon integration against 5(u ■ x), we obtain 

5Aare = ' T Pb(^) -^q^T. (A.20) 

J Tl 

Here, Pb is the Minkowski momentum four-vector of the bare particle, which for later con- 
venience we write in the form 

Pb = Mb-u, (A.21) 

where 

Mb = -Mbg (A.22) 

is the diagonal bare mass tensor, with 

-Mb(||w,||)= / ^ /^(||x||)^(u-x)d^x (A.23) 

■^^^ ^Jl-\\n^■x\\^ 

the gyrational bare mass of the particle. 

Turning next to the action of field-particle coupling, we rewrite it as 

Aoupi = / ^c(f^E, q, ", ") dr , (A.24) 
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where 

Lc(nE, q, u, u) = /" (u-n-x)- A(q + x)/e( ||x|| )5{u ■ x) d^x (A.25) 

is the Lagrangian for the couphng of the particle to the electromagnetic fields. Since fig and 
u occur in L^ only in the combination fig + u A u = f2, we can treat L^. as a function of 
n, q, u. Thus, abusing notation a little bit, we write Lc{0,, q, u) and obtain 

^Aoupi = f (^ Tr(^ ■ <50) + ^ ■ <5q + ^ ■ <5u - L,u ■ /q)dr (A.26) 

under independent variation of the world-line. As for the variation of Lc w.r.t. u, we have 

6Lr 



6 



-= I A(q + x)/e(||x||)5(u-x)d^x 

J-.q 

+ I {u-n-x) ■A(q + x)/e(||x||)9„(5(u-x)d^x. (A.27) 

J ^Q 



Using once again (|A.16|) , then integrating by parts, carrying out the differentiations, and 
noticing once again that terms proportional to u ■ x in the integrand vanish upon integration 
against (5(u ■ x), we obtain 

^= / A(q + x)5(u.x)/,(||x||)d^x 

+ u /" {u-n-x)- A(q + x)/e( ||x|| ) 5{u ■ x) d^x 

J'B.q 

+ I xu- A(q + x)/e(||x||)5(u-x)d^x 

J=.q 

+ I x{u-n-x)-{u- VgA(q + x))/e( ||x|| ) 5{u ■ x) d^x. (A.28) 

Using the four-orthogonality of u and 5u, we thus have 
^-^.5u= I 5u-A{q + x) 5{u ■ x)/e( ||x|| ) d^x 

+ / (5u ■ xu ■ A(q + x)/c( ||x|| ) 5(u- x) d'^x 

J'Eq 

+ I 5u-x{u-n-x)-{u- VgA(q + x))/e( ||x|| ) 5{u ■ x) d^x. (A.29) 
As for the variation of Lc w.r.t. q, we have 



^^^ 6q = Jju-n.x). {^^^^§^ ■ 5q^ /e( l|x|| ) 6{u ■ x) d^x, (A.30) 



6q 

50 



while for the variation of L^. w.r.t. fi we simply have 

- Tr(-^ ■ 5Vl) =- I i6n ■ x) ■ A(q + x) /e( ||x|| ) 6{u ■ x) d^x. (A.31) 



We use ( |A.7] ) for the variation of the total gyration tensor f2 along the independently per- 
turbed particle world-line, and find 



SA, 



coupl 



/T2 J I A r 

— I -(g+u ®u) ■— / xu- A(q + x)/e(||x||)(5(u-x)d^x 

+ / [A(q + x) - (u • A(q + x))n, ■ x] /e( ||x|| ) 6{u ■ x) d^x 

+ f [x+(u-x)u](u-0-x)-A(q + x)/e(||x||)(5'(u-x)d^xUqdr 
Vg® A(q + x) • [u-rj-x]/c(||x||)5(u-x)d^x5qdr. (A.32) 



We next simplify the previous equation with the help of the following useful identity. Per- 
forming an integration by parts, one can easily check that 

f [u - O • x] • Vg ® A(q + x)/e( ||x|| ) S{u ■ x) d^x 
= [ A(q + x)(Trf2)/e(||x||)5(u-x)d^x 

- f A(q + x)[u-x-x-fi-x]/^(||x||)/||x|| 5(u-x)d^x 

- f A(q + x)u- [u-0-x]/e(||x||)5'(u-x)d^x. (A.33) 

By the anti-symmetry of O we have Tr O = 0, so that the first integral of the r.h.s. of ( [A.33| ) 
vanishes. Again by the anti-symmetry of ft, we also have that x ■ f2 ■ x = 0, and furthermore 
we have u- x6{u- x) = 0. Therefore, the second integral of the r.h.s. of ( [A. 3^ ) vanishes, too. 
As for the last integral of the r.h.s. of ( |A.33|) , we use that f2 = fig + f^Fw with Op^ = u A u, 
and moreover that u ■ fig = 0, to rewrite it into 

- I A(q + x)u- [u-ri-x]/e(||x||)(5'(u-x)d^x 

f A(q + x)[l + u-x]/e(||x||)5'(u-x)d^x. (A.34) 

JSq 

Finally, the r.h.s. of ( |A.34D is integrated by parts once, and we arrive at the final identity 
f [u - f^ ■ x] ■ Vg ® A(q + x)/e( ||x|| ) 6{u ■ x) d^x 

= ^ [ A(q + x)/e(||x||)<5(u.x)d^x. (A.35) 



51 



This identity is now used to anti-symmetrize the last integral in ( |A.32|) . Noticing next that 



Vg A • V — V ■ Vg ® A = Vg A A ■ V for an arbitrary four-vector v, and recalling that 

coupi takes the form 



Vg A A = F, the variation ( [702D of A 



SA. 



coupl 



T2 



p^ - f Sq dr 



where 



f=f F(q + x).(u-f^-x)/e(||x||)(5(u-x)d^x 
is the Abraham-Lorentz type Minkowski force, and 



(A.36) 



(A.37) 



Pn 



X • f^E ■ F(q + x) ■ u /c( ||x|| ) 6{u ■ x) d^x 



(A.38) 



is Nodvik's Minkowski momentum for the electromagnetic spin-orbit coupling. For later 
convenience, we notice that the charged particle's co-moving electromagnetic field, viz. its 
Coulomb field and its instantaneous gyro- magnetic field (both defined w.r.t. the Fermi- 
Walker frame) do not contribute to p^, so that F in ( |A.38|) can be replaced by the reduced 
field tensor F', which obtains from F by subtracting the particle's co-moving electromagnetic 
field. Furthermore, making use of the identity u ■ fig = 0, and using the fact that terms 
proportional to x- u in the integrand vanish upon integration against 6{u-x), we can rewrite 
the four- vector p^ in the more appealing format 



Pn 



M. 



u 



where we introduced the symmetric Nodvik tensor mass of spin-orbit coupling. 



M, 



X ® X, [F'(q + x), n^] /e ( ||x|| ) (5(u ■ x) d^x . 



(A.39) 



(A.40) 



The Euler-Lagrange equation for 6{Aha.Tc + ^coupi) = under the variation of q(r) for 
given F(x) and fl-Ei^) can now be read off from ( |A.36| ) and ( [A.20| ) as 



dp 

dr 



where 



M u 



(A.41) 



(A.42) 



is the four- vector of the Minkowski momentum, with M = Mb + M^. Equation (|A.41| ) is the 
Minkowski equation of motion with the expected Abraham-Lorentz type Minkowski force, 
but with a somewhat unusual anisotropic tensor mass M. 
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A. 1.4 Euler-Lagrange equations for the gyrograph 



Using (|A.8| ) and integrating by parts, observing that the variations of the Euler variables 



vanish at the boundaries, for the first variation of the action Aha.ve of the bare particle we 
find 

SA^.re = \ r Tr ((Sb - [Sb M]_)- so) dr, (A.43) 

and for the first variation of the coupling action we get 

5Aoupi = I P Tr ( (Sf - [Sf , fl] _) • 56) dr (A.44) 

where we introduced the tensor of electromagnetic Minkowski spin (about q) of the charged 
particle, 

Sf = /" X A A(q + x)/e(||x||)5(u-x)d^x. ( p39D 

We remark that Sf is invariant under a gauge transformation A ^ A + Vg\E'. Defining now 
the total spin of the particle as 

S = Sf + Sb (A.45) 

and introducing the abbreviation 



def 



N = ^S + [0,S]_, (A.46) 

dr 



the vanishing of the first variation of the action functional with respect to all anti-symmetric 
60 that are of space-space type with respect to u now simply reads 

0= ['tt{N-60) dr. (A.47) 



At this point we recall that, given any tetrad, any anti-symmetric tensor A can be 
uniquely decomposed into a sum of its space-space and time-space parts. With u as the 
timelike unit vector of the tetrad, we have 

A = A^ + a" , (A.48) 

where 

A"*" = A + [u (g) u , A] is the space-space part of A, and , , 

a" = u a (A ■ u) its time-space part. 



Armed with the above decomposition, we conclude that equation ( |A.47| ) can only be 



satisfied if the space-space part w.r.t. u of the anti-symmetric tensor N vanishes, explicitly 
N^ = 0, (A.50) 
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which is the Euler-Lagrange equation for the rotational variations. 

To extract the evolution equation for flj, from (|A.5CI| ), we also decompose the anti- 



symmetric spin tensor S into its space-space and time-space parts with respect to u, 

S = S^ + Sll. (A.51) 

Inserting the above decomposition of S into the definition ( |A.46| ) of N simply gives 



N = ^S^+[0,S^1 +^Sll+[0,Sin . (A.52) 

(XT '- -' " dr '- -' ~ 

With the help of the helicity four- vector Jf = S • u we can write 



S^ = S - u A H 

S" = u A H 

By straightforward calculation we furthermore find that 



(A.53) 



d7 



s" + [n, s"] _ = ua(h+ n^-Hj 



and that 



{is-+[n,s-]_yu = o. 



(A.54) 



(A.55) 



Recalling next that S^ ■ u = (i.e. Sb" = 0, so that Sb = Sb"*") we find that S ■ u = Sf ■ u, 



VIZ. 



/ X u- A(q + x)/e( ||x|| ) (5(u- x) d'^x 



H= xu- A(q + x)/e(||x||)5(u-x)d^x, (A.56) 

from where we read off that H ■ u = 0, i.e. H is itself spacelike for all r. This implies that 

o 

Jf + fig ■ Jf is spacelike, and we conclude that 

^Sll+[fl,Sl'l =Nl', (A.57) 

(XT '- -' 

^S^+rn,S^l =N^. (A.58) 

dr 

Hence, the Euler-Lagrange equation ( [A. 501 ) reads 



— S^+ rO,S^l =0. (A.59) 

dr '- -'^ 

Remark: The conservation law for IIS 11 follows immediately from ( |A.59| ). 

It remains to show that ( |A.59|) is identical to ( [7.6|) . We recall that S = Sb + Sf, and 



notice that decomposition into space-space and time-space parts commutes with addition, 
so that S = Sb -l-Sf . We also recall that Sb = Sb. We furthermore notice that fig oc Sb, 
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which imphes that [OE,Sb]__ = 0, so that [f2, Sb]_ = [f2Fw,Sb]_. Using all these pieces of 
information we see that (|A.59|) states that 



^Sb + [flpw, Sb] _ = -^Sf^ - [n, Sf^] _ . (A.60) 



The l.h.s. of ( |A.60[ ) is manifestly identical to the l.h.s. of (|7.6| ). To see that the r.h.s. of 



( |A.6CI| ) is identical to the r.h.s. of ( [7.6| ), one can rewrite the calculation in Appendix C of ||T5[ 
into our notation, using the Nodvik Maxwell-Lorentz equations (|5.1|) , (|5.2|) , (|5.3|) , and a few 
integrations by parts. The details are not repeated here. 

A. 2 Point mass and point charge limits of massive LED 

In the renormalization flow study of our massive LED with the choice (|10.1|) for the mass 
and charge densities, the renormalized particle reaches a finite smallest size = 1.5 -Re in 
the 'renormalized purely electromagnetic' limit of vanishing bare rest mass mb -^ 0+. This 
renormalization program for massive LED aside, it is of interest to investigate which types 
of theories emerge in the point limits for mass and charge, either taken separately or jointly. 

A. 2.1 The point mass limit 

In the point mass limit 

UM) — rn,-^ 6{\\x\\) '^ /r°^-( ||x|| ) , (A.61) 

but with /c fixed and 'regular' (not a point charge density), and with the condition that 



||ive|| < C, our massive LED formally reduces to the massive LED of Nodvik [T^. In this 
limit, the square root terms in ( [4.2| ) and ( |6.5| ) become just factors of unity on the support 
of f^'^^^^^, so that the bare particle's gyrational mass ( ^4.2| ) degenerates identically into a 
constant, 

A^bdlw^lD^mb, (A.62) 

and our bare action reduces simply to the action of a free mechanical point particle with 
bare rest mass mb, 

£bare(x) d^x -^ -m^ H dr = AZT"" , (A.63) 

Jti 



employed by Nodvik [|1^]. The field and coupling actions remain the same, however. Ac- 
cordingly, the electromagnetic part of the two models is identical, i.e. the Maxwell-Lorentz 
equations in Nodvik's massive LED are given by (^.1|), (|5.2|), and (|5.3|). Moreover, the 



d_ 



world-line equation in Nodvik's massive LED can still be written as 

-p{r)= [ F(x)-C/(x;r)/e(||x-q(r)||)<5(u.(x-q(r)))d^x, (A.64) 
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but now one has 

p{t) = MN°^"''^(r) ■ u{t) (A.65) 

with 

MNodvik(^) ^ ^^g _ /■ [x^x, [F(q + x),i7E]J /e(||x||)(5(u.x)d^x, (A.66) 

which is a shghtly simpler mass tensor than M. On the other hand, the point mass hmit 
ehminates the moment of inertia of the bare particle while the norm of fig remains bounded 
(by hypothesis) so that Sb — > identically; whence, the gyrograph equation (|7.6| ) degenerates 
into the statement t = 0, i.e., upon using ( |7.7| ), 







/ (x-q)A(F(x)-U(x,r))Ve(||x-q||)<5(u.(x-q))d^x. (A.67) 

jRl.3 



These are the equations of the massive LED of Nodvik [|I5| . We will now explain that this 
is a qualitatively quite different dynamical theory from our massive LED. 

First of all, since one loses a proper-time derivative for fig by passing from ( |7.6|) to ( |A.67|) , 
the point mass limit is a singular limit from the dynamical perspective of the Cauchy problem 
for rig. This fact in itself already means that it will no longer be possible to arbitrarily choose 
initial data for fig from the set of admissible Cauchy data for (|7.6| ). Indeed, choose (which 
we always may) a Lorentz frame in which initially (i.e. at t = = r) we have qg = and 
uq = Cq, and pick the fields Fq on uq ■ x = at random from our set of admissible Cauchy 
data for Fq. Then ( |A.67 ) at r = has only a nontrivial space part, which in three-space 



notation reads 

t^ + cj^oxa^ = 0, (A.68) 

where uj^q is the space part of the dual four- vector for JIeQ) ^^d where 

t^= [ xx£o{x)U\x\)d'x (A.69) 

is the initial electric torque on the particle and 

^0^^=' / xx-Bo{x)U\x\)d'x (A.70) 



its initial 'magnetic spin.' Clearly, any attempt to also prescribe Cauchy data for the angular 
velocity tensor fl^o, chosen at random from our set of admissible data, will almost surely 
violate ( |A:6^ ). 

More specifically, ( |A.6^ ), with ( |A.69|) and ( |A.70|) , imposes not just on the set of admissible 
initial data for fig but also on the fields Fq, in the sense that for most choices of Fq no choice 
for J7eo will satisfy (|A.68|) . Thus, one extracts directly from (|A.68|) , ( |A.69| ), (|A.70|) that 



the originally admissible field data are now further restricted as follows: either Bq satisfies 
cr^ = 0, and in this case also £q has to satisfy Iq = 0; or Bq satisfies ctq ^ 0, in which 
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case Eq has to satisfy t^ 



0. Only if these alternate constraints on the initial fields are 



satisfied is it possible to satisfy (|A.68| ) by a compatible class of initial data f^EO (viz-, ^eo) 
Namely, if cTq = and t^ = 0, then the admissible data for cjpn are the same as for our 
model. 



'EO 



If <T^ ^ and t^ 



E 



0, however, then ujpn has to be of the form 



(jj 



EO 



rB\2 



'-0 



^? 



t^ X an + acr 



(A.71) 



which leaves only one degree of freedom to choose from instead of three, i.e. the coefficient 
a in front of a^ is at our disposal. This clarifies to which extent initial data for Nodvik's 
model have to be compatible. 

Compatibility of a set of initial data in the sense just explained does not, however, imply 
that all is well then, for not only has one lost a proper-time derivative on Q,^ by passing 
from ( [7.6|) to ( A. 67 ), a proper-time derivative of Oe does not show any more at all in ( A. 67 ). 
Nodvik's ( |A.671 ) is therefore not an orderly evolution equation for f^E, but rather a constraint 
equation for the dynamical field variables F and the dynamical particle variables q and u, 
which in case of actually being satisfied by these variables for some time interval after the 
initial instant, still leaves a certain amount of freedom for choosing the value of the now 
purely kinematical variable f^E- Therefore, to the extent that compatible initial data for 
F, q, and u, with a compatible initial choice of f^E, can be continued at all according to 
Nodvik's equations, such a continuation is presumably not unique. 

Remark: We shall see that a similar assessment holds for the semi-relativistic, purely 
electromagnetic theory of Abraham, which is discussed in Appendix A. 3. 



A. 2. 2 Remarks on the spinless point particle limit 

Ever since Dirac's surmise [0] that, as long as we are willing to let the bare mass take 
the value — oo to compensate for the +oo electromagnetic self-energies, we may view the 
Abraham-Lorentz-Dirac (ALD) equation 



00 — .t° '^ ^ / ° ° \ °°° 

mq = -eF°^' ■ q + - — (g + q (g) q) ■ q 
o c 



(A.72) 



as the exact equation for the radiation reaction-corrected motion of a spinless point particle 
with empirical electron mass m and empirical electron charge — e in external electromagnetic 
fields F™*, there has been a strong interest in making sense out of the peculiar third-order 
proper-time derivative of q that figures in ( A.72|) , e.g. |27, [ 



62, 33, 64 



|. A rigorous 

clarification of Dirac's surmise is, however, not so easy to come by. Notable rigorous attempts 
to construct a consistent renormalized electrodynamics with spinless point charges are |Jl9| 
and [Q; see also [^ and |^, ^ for work on the classical Pauli-Fierz model, a semi- 



relativistic dipole approximation to the problem, which also serves as point of departure for 



setting up a so-called non-relativistic QED ||53 |. 

Naturally one may wonder whether ( [A.72[ ) emerges in some spinless point particle limit 
of our massive LED, or to begin simpler, in the point particle limit of a spinless version of our 
massive LED. This massive LED without spin obtains from our model by simply discarding 
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(7.6) and setting ilg = in the remaining equations. Thus, the relativistic Maxwell-Lorentz 



equations are once again 

Vg • '^F = (A.73) 

for the Hodge dual of F, and 

Vg ■ F(x) = 47r /" " l/(x; r) /e( ||x - q(r)|| ) 5(uir) ■ (x - q(r))) dr (A.74) 

J —oo 

for F itself, but now with 

l/(x; r) = u{r) - {uir) A u(r)) ■ (x - q(r)) . (A.75) 

The world-line equation still reads 

-p= [ F(x)-U(x;r)/e(||x-q(r)||)5(u(r)-(x-q(r)))d4x, (A.76) 

but now with 

p{T)=m^u{T), (A.77) 



d 



and with U given in ( [A.75[ ). Notice that the structure of /m does not figure in this spinless 



LED, which is therefore identical to the spinless version of Nodvik's LED, due to Spohn [^ 



in which f^ is a point mass distribution from the outset. The question then is: What 
happens if we let /e tend to a delta function concentrated at the same point? The first step 
to the answer is the rigorous control of this spinless massive LED which, however, is not 
very developed yet. 

A. 3 Semi-relativistic LED 

In this Appendix we give a brief summary of semi-relativistic LED. This material can also 
be used as a primer for the relativistic theory. We begin with the massive model with spin. 
Its limiting cases of the more familiar Abraham model and the popular 'model without spin' 
are discussed in subsequent subsections. 

A. 3.1 Massive model with spin 

As in the relativistic theory, the electric field £{x, t) G M'^ and the magnetic field B{x, t) G M? 
at the point x G M'^ at time t G M satisfy the classical Maxwell-Lorentz equations, 

V-B{x,t) = Q, (A.78) 

-^jB{x,t) + V^£{x,t) = Q, (A.79) 

V ■£{x,t)=A-Kp{x,t), (A.80) 

Id 1 

— ^£{x, t) + V X Bix, t) = Anpix, t)-v(x, t) , (A.81) 

cat c 
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Although not manifestly Lorentz co- variant, the Maxwell-Lorentz equations ( |A.78| ), ( A.79|) , 
( |A.8CI| ), (|A.81|) appear to have a relativistic format already, but whether in fact they do 
satisfy relativistic Lorentz transformation laws or not depends on the model for the charge 
distribution. In the classical works of Abraham |jl], ^ before the dawn of special relativity 0, 
^, H, H, the electrical charge and current densities formally combine as p{x, t) (c, v{x, t)) but 
do not transform like a relativistic four-current density vector. 

In the semi- relativistic formulation of Abraham, the charge of a particle (— e for a model 
of the electron) is distributed around the instantaneous location q{t) G M.^ of the particle by a 
charge density /e with S0{3) symmetry and compact support, satisfying J^^ /e(|a;|)d^x = — e. 
This charge density is rigidly carried along by the particle with linear velocity q{t) and 
rotating rigidly with angular velocity uj{t). Thus, with a single charged particle source in 
( lOOl ), ( |Ol| ), we have 



p{x,t) = U\x~qm, 



(A.82) 



and 



v{x, t) = q{t) + uj{t) X [x- q{t)) 



(A.83) 



The evolution equations for the dynamical variables of the particle, i.e. position q{t), 
linear velocity q{t), and angular velocity ci;(t), are likewise not Lorentz covariant and given by 
Newton's and Euler's evolution equations equipped with the Abraham-Lorentz expressions 
for the volume-averaged Lorentz force and torque |^ , respectively, that act on the charged 
particle. Newton's equation of motion here reads 



dt 



S(x,t) + -v(x,t) X B(x,t) 
c 



M\x-qmd'x., 



where 



(A.84) 



Pb = m^q 



(A.85) 



is the Newtonian linear particle momentum, with r/ib the bare mass ('material mass' in ||13|| ). 
Euler's equation of motion reads 



dsb 

dt 






where 



£{x,t) + -v(x,t) X B{x,t) 
c 



U\x-qmd'x., 



(A.86) 



Sb = h^ 



(A.87) 



is the classical particle spin associated with the bare moment of inertia Jb. 

With the bare inertias ?7ib 7^ and Jb 7^ 0, the semi-relativistic equations listed above pose 
a Cauchy problem for the following initial data at time 1 = 1^. for the mechanical variables 
of the particles, the data are g(to), qit^) (equivalently: Pb(^o)); and ijjito) (equivalently: 
Sb(to)); for the fields, the data are B{x, to) satisfying ( |A.78| ) at t = to; and S{x, tg) satisfying 
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( A. sop at t = to- Notice that (|A.78|) and ( |A.80| ) are merely initial constraints, which then 
remain satisfied by the fields B{x,t) and E{x,t) in the ensuing evolution. For (|A.78|) this is 
seen by taking the divergence of ( |A.79| ). For ( [A.80| ) this is seen by taking the divergence of 
( |A.81| ) and the time- derivative of ( |A.80|) , then using the continuity equation for the charge, 



dtpix, t) + V ■ {p{x, t)v{x, t)) = 



(Ai 



which is demanded by the Maxwell-Lorentz equations and indeed satisfied for p{x,t) given 
by ( pm and v{x,t) by (|X:83|) ; cf. ^. 

The semi-relativistic massive LED with spin has the following conservation laws 
Charge, 



p{x,t) d^x . 



energy. 



2 r/ib 



linear momentum, 

V = -^ [ SxBd?x + p^, 

and angular momentum, 

C = / X X (S X B) d^x + qr X pjj + Sb 

4vrc Jks 



(A.89) 



(A.90) 



(A.9i: 



(A.92) 



are all conserved during the evolution. 

More detailed studies of the coupled dynamics of the massive semi-relativistic LED with 
spin have only recently begun. 

A. 3. 2 Massive model 'without spin' 

A mathematically regular limit which, unfortunately, is not quite physical but brings about 
some considerable simplification of the dynamical equations is to send Jb ^ oo while rri]^ 
is fixed at some finite value, in such a way that Sb converges to some finite vector- valued 
function of t, which of course implies that a; ^ 0. The semi-relativistic Maxwell-Lorentz 
equations are again given by (|A.78|) , ( |A.8CI| ), (|A.79|) , ( |A.81[ ), and the charge density is still 
( |A.82| ), but the current density vector now simplifies to 



j{x,t) = p{x,t)q{t). 
Moreover, the Newtonian equation of motion ([A.102|) reduces to 



dPb 
dt 



£{x,t) + -q{t) xB{x,t) 
c 



U\x-qmd'x., 



(A.93) 



(A.94) 
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with p still given by ( |A.85|) , and the Eulerian equation ( |A.103| ) reduces to 



dsb 
dt 



I [x- qit)) X 

JR3 



£{x,t) + -q{t) xB{x,t) 
c 



U\x-qmd'x. 



(A.95) 



In this limit, the Maxwell-Lorentz equations ( |A.78| ), ( |A.80| ), ( [A.79| ), ( |A.81| ), equipped with 
the charge density ( A.82|) and current density vector ( |A.93| ), and the Newtonian equation 
of motion ( |A.94| ), form a non-linearly coupled, closed system of equations, with initial data 
at t = to given by: q(to) and q'(to) (equivalently: Pb(^o)) for the particle; for the fields, 
the data are B{x,tQ) satisfying ( |A.78D at t = to, and £^(a;,to) satisfying ( |A.8CI| ) at t = to- 
Spin in turn adjusts passively to this combined field-particle dynamics according to a linear 
equation ( |A.95| ), which has to be supplemented by the initial data Sb(to) (now of course 
no longer equivalent to prescribing a; (to)). Consequently, in this limit spin drops from the 
conservation laws — with the exception of the law for the angular momentum. Thus, the 
conserved energy W, linear momentum V, and angular momentum C are now given by 




'x + - 



llPbl 



2 m\j 



(A.96) 



Anc 



and 



Anc 



£ X B d^x + Pb , 



X X {£ X B) d^x + q X p^ + S),. 



(A.97) 



(A.98) 



Notice that the laws for energy and linear momentum do not feature Sb? but the conservation 
law for the angular momentum does. Furthermore, the continuity equation simplifies to 



dtpix,t) + V-{pix,t)qit)) =0 



(A.99) 



but this does not affect the law of charge conservation, which is still given by ( |A.89| ). 

Since, in the limit /b — > oo, spin is of no dynamical relevance, apart from the law of 
angular conservation, ( |A.95| ) is frequently omitted from the system of equations. The result- 
ing subset of equations defines the semi-relativistic massive LED 'without spin.' This model 
without spin is still invariant under rotations, and by E. Noether's theorems there exists a 
conserved quantity associated with this invariance; only this quantity does not qualify as the 
physical angular momentum in the conventional sense. The observation that discarding the 
gyroscopic equation ( |A.95| ) violates the conventional law of angular momentum conservation 
( |A:98| ) was made in l7|. 

The semi-relativistic massive LED without spin is the most thoroughly understood ver- 
sion of LED. In particular, we mention the recent rigorous treatments in |7^, |73|, 0, |75|, [76 
where for technical reasons the Einsteinian linear momentum 



Pb 



m^q 



|qr|2/c^ 



(A.lOO) 
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with r/ib 7^ is used instead of the Newtonian (|A.85|) (see also |77, 78 1 for work on a simpler 
scalar field-particle theory). In this case, of course also the Newtonian kinetic energy in 
( |A.9CI| ) is to be replaced by its Einsteinian counterpart 



1 It) |2 

-*- IPbl 

2 rrih 



m^c 



1 + 



iPbl 



mtc 



2„2 



(A.lOl) 



In |7^, |73] , the global existence and uniqueness for the Cauchy problem of a particle without 



spin and m\^ ^ was proven. Moreover, in [^ it was shown that the motion is stable if 
mb > and unstable if rrib < 0. The papers |7^, |7^, |7^, |76|] address the long time asymptotics 
of the slow motion of a spinless particle with mb > in slowly varying external fields. Using 
center manifold theory, these authors derive effective equations of motion of second order 
in the time derivative. Originally these effective second-order equations were obtained in 
Landau-Lifshitz [|l^ by applying a heuristic closure argument to the (in)famous third-order 
term in the Abraham-Lorentz-Dirac (ALD) equation ( |A.72D , which works equally well for 
the relativistic and for the semi-relativistic version. As far as we know, the rigorous work of 
Spohn and collaborators gives the first clean derivation of the Landau-Lifshitz closure for 
the ALD equation, so far in its semi-relativistic setting. 



A. 3. 3 The Abraham model with spin 

Abraham's [jl], 0] purely electromagnetic model obtains from the massive semi-relativistic 
LED with spin in the formal double limit rrib — > and /b — > together with the assumption 
that the velocity q and angular velocity oj of the particle remain bounded. The Maxwell- 
Lorentz equations are not affected by this limit and still given by (|A.78|) , (|A.79|) , (|A.8CI|) , 
( |A.81| ) together with ( [A.82| ) and ( |A.83| ). However, the bare Newtonian momentum Pb ~^ 
and bare Eulerian spin Sb ^ identically in this limit. As a consequence, the evolution 
equations ( |A.84D and ( |A.8(j| ) degenerate into 



and 











Six.t) + -vix,t) xB(x,t) 
c 



U\x-qmd'x 



{x - q{t)) 



£(x,t) + -v(x,t) X B(x,t) 
c 



U\x-qmd'x. 



(A.102) 



(A.103) 



The conservation laws for energy, linear and angular momentum simplify accordingly, be- 
cause Pb ^^^ *b vanish quadratically fast with m-b —>■ and Jb — > while q and u: remain 
bounded (by hypothesis!). Thus, the conserved energy W, linear momentum P, and angular 
momentum C are now given by 






X . 



(A. 104) 



Aire 



£xB& 



X . 



(A.105) 
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and 

/: = ^ I XX (E xB)d^x. (A.106) 



The continuity equation ( |A.88| ) and its offspring, the law for the charge ( |A.89| ), still hold as 
stated. 

This purely electromagnetic Abraham model with spin is not well-studied at all. In fact, 
it is hardly mentioned in textbooks and monographs on classical electromagnetism, with the 
notable exception of the book by Miller [Q . 

Subjecting equations ( |A.102| ) and (|A.103|) to closer scrutiny reveals that Abraham's 



purely electromagnetic semi-relativistic LED with spin is qualitatively a very different dy- 
namical theory from the massive semi-relativistic LED with spin. In fact, since the coeffi- 
cients in front of the time derivative of both q and a; vanish, these equations constitute a 
singular limit of the regular evolution equations (|A.84 ) and (|A.86 ), a limit in which Cauchy 



data for the dynamical particle variables q and lj can no longer be prescribed freely; also 
the field data are severely restricted. Only the initial particle position can be prescribed 
arbitrarily in Abraham's model. The situation is similar as in Nodvik's relativistic model, 
only worse. 

More specifically, choose the initial time to be t = and the initial particle position to 
be q{0) = 0, which can always be achieved by at most a time and a space translation. As 
before, denote the initial fields by Soi^) and Bo{x). Then ( [A.102| ) demands that at t = 0, 
we have 

c{£o) + go X (^o) -{{x(^Bo)- (tr {x®Bo))\)-u;o = 0, (A.107) 

while (|A.103 ) demands (as for Nodvik's model) that at t = we have 



c{xx£o) +ujoX (xx-Bo) = 0. (A.108) 

Here, the averaging is defined by {g) = {—e)^^J^3g{x)fe{\x\)d^x, where we assume with 
Abraham that (— e)"^/c(|a3|) ^ 0; furthermore, Cg> now means the three-tensor product. One 
readily deduces from ( |A.107| ) and ( |A.108| ) that only the following subset of the originally 



admissible field data is now admitted. Namely, if {xx ■ Bq) = 0, then also {x x Sq) = 
must hold; if (a; a; ■ Bq) ^ 0, then {x x £q) ■ {xx ■ Bq) = must hold; if {B^) = 0, then {Eq) 
must be in the orthogonal complement of the left kernel space of {x ® Bq) — (tr {x ^ Bo))\, 
unless this tensor is the zero tensor, in which case (So) must vanish also; if (Bo) 7^ is in 
the left kernel space of {x (g) Bq) — (tr {x <S) Bo))\, then also (£0) ■ (^0) = must hold. 

Remark: The above set of conditions on the initial fields imply in particular that the 
Abraham equations admit no solution for the following setup: q^ = 0, Bq = 0, and Eq = 
£hom_|_£ Coulomb^ j ^ ^j^^ linear superposition of the charge's radial symmetric Coulomb field 
with a homogeneous external field. In this case no choice of ^q ^^^ '^0 will rescue ( |A.107|) ; 



in short, these perfectly physical initial conditions are inconsistent for the Abraham model. 

We continue with the discussion of consistent initial data. Given now consistent data for 

the fields in accordance with the restrictions just explained, the initial data for qg ^^^ ^0 
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are arbitrary only if (Bq) = and {x x ■ Bq) = 0, respectively; otherwise they must be of 
the following form. If {xx ■ Bq) 7^ 0, then 

{x X So) X {xx ■ Bo) 



Uo 



+ -f{xx-Bo) 



(A.109) 



\{xx-Bo)\' 

which leaves only one degree of freedom to choose from instead of three, viz. the coefficient 
7 is at our disposal. If {Bo) 7^ 0, then 

1 



Qo 



i£o) + {{x (g) Bo) - (tr {x (g) Bo))\) ■ cuo X (^o) + f3{Bo) 



(A.llO) 



\{Bo)\' 

which again leaves only one degree of freedom to choose from instead of three, here the 
coefficient f3 is at our disposal. The extent to which such consistently restricted initial data 
now launch a unique regular solution of the system of Abraham equations ([A.102| ), ( |A.103|) , 
( 1^:751 ), (^7^), ( |OIl| ), and (|X:HT[ ), with (^[^ and (^[S^), is currently not known. 



A. 3. 4 The Abraham model without spin 



The purely electromagnetic Abraham model without spin [|79| , |80| obtains from the Abraham 
model with spin by simply discarding ( |A.103| ) and setting a; = identically in ( A.83|) , so 
that ( [A.102|) goes over into 







e{x,t) + -q{t) xB{x,t) 



fMx-qmd'x., 



(A.lll) 



while the Maxwell-Lorentz equations for this model are the same as for the massive model 
without spin, i.e. again given by ( |A.78| ), (|A.80|) , (|A.79| ), ( |A.81|) , with the charge density 



( |A.82| ) and with the current density vector ( [A. 93] ). The continuity equation for the charge 
is the simpler (|A.99|) , while the law of charge conservation is of course just ( |A.82| ). The 
conserved energy, linear and angular momentum are given again by the purely electromag- 
netic expressions ( |A.104D , ( |A.105| ), and ( [A.106| ) that apply to the Abraham model with spin. 
The Abraham model without spin can also be obtained as the singular limit mb ^ of the 
massive semi-relativistic model without spin under the hypothesis that q remains bounded. 

For the Abraham model without spin a similar assessment holds as for the Abraham 
model with spin. Again choosing the initial time to be t = and the initial particle position to 
he qo = (after at most a time and a space translation), the restrictions for the field data are 
now determined by ( JA.IOTD with ujo = (and ignoring the discarded ( |A.108D ). We conclude 
that, if {Bo) = then also {£0) = must hold, and if {Bo) ^ 0, then (£0) ■ (^^o) = must 
hold. Although these are milder restrictions than those for the Abraham model with spin, 
these restrictions on the field data in the Abraham model without spin still rule out various 
perfectly physical initial data, for instance those remarked on in the previous subsection, i.e. 
q^ = 0,Bo = 0, and £0 = fhom ^£Coulomb_ 

As for the initial velocity of the particle q^ in given consistently restricted initial fields, 
the following rules hold. The initial velocity q^ is at our disposal only if {Bo) = 0, and if 
{Bo) 7^ it is of the form 

{So) X {Bo) 



Qo 



\{Bo) 



+ mo) , 



(A.112) 
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which, once again, leaves us only the choice of p. 

Whether such a restricted, consistent set of initial data can now be continued into a 
unique regular solution of (^TTTl|), PTTSP , piSOP , (^T|), and d^ISTD , with (^^) and 
( |A.93| ), is currently not rigorously known. As a matter of fact, there are good reasons to 
doubt that the evolution problem is well-defined even given consistent initial data. For notice 
that ( |A.lll ) states that at 



{sm + -qit)x{Bm 



0. 



for all t, where the averages here are now defined by 



{gm 



—e 



g{x,t)f,{\x-q{t)\)d'x. 



(A.113) 



(A.114) 



Hence, we can repeat our discussion of the conditions on the initial data essentially verbatim 
to conclude that: 



If at any time t we have {B){t) = 0, then ( |A.113| ) has no solution q{t) at that instant 
t unless also {S)(t) = 0; 



If {13){t) ^ 0, then ( [ATT^ ) has no solution q{t) at t unless also {£){t) ■ {B){t) = 
If {B){t) ^ and {€){t) ■ {B){t) = 0, then (|A.113| ) states that q{t) is of the form 

{sm X {Bm 



m 



mitw 



+m{Bm 



(A.115) 



with an a priori undetermined coefficient l3(t). 



It does not seem very likely that this system of equations, coupled to the Maxwell-Lorentz 
equations, evolves a unique regular solution from a consistent set of initial data. 

While there is hardly any literature on the Abraham model with spin, the Abraham model 
without spin is discussed in most major texts since Abraham 
and also in the more leisurely text of Peierls 



e.g. ig, n, n, grm, |T 



] . However, none of these texts makes the el- 
ementary observation that the initial value problem for the purely electromagnetic Abraham 
model(s) is singular in the sense explained above, nor that Abraham's 'evolution' equation 
( |A.113| ) preserves this singular structure in time; some of the main conclusions drawn in 



[^, ^ ^7\, q1|, ^ are even in open conflict with the singular nature of Abraham's evolu- 
tion equations. A re-assessment of these traditional treatments of the Abraham model seems 
to be called for. 
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